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ABSTRACT 


The Rocketdyne Safety Algorithm (RSA) has been developed to the point of use on the 
TTBE at MSFC on Task 4 of LeRC contract NAS3-25884. This document contains a 
description of the work performed, the results of the nominal test cases and a table of 
the resulting cutoff threshold levels, a plot of the RSA value vs. time for each nominal 
case, the results of the major anomaly test cases and a table of the resulting cutoff times, 
a plot of the RSA value vs. time for each anomaly case, a logic flow description of the 
algorithm, the algorithm code, and a development plan for the future efforts. 
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1.0 SUMMARY AND CONCLUSION 


The Rocketdyne Safety Algorithm (RSA) has been developed to the 
point of use in a test environment. The objective of this program is 
to develop the algorithm to a form which can be tested against 
previous engine test data sets for validation and coded for use with 
actual engine test data at MSFC. The real time influence coefficient 
engine model has been improved to the level that it works 
throughout the SSME sequence from start to cutoff. As a result, the 
RSA is functional throughout a complete test, during transients as 
well as steady state operating periods. The algorithm has been coded 
in Ada and verified to work on the SAFD hardware. However, due to 
timing constraints in the SAFD system, it will not be able to run 
concurrently with the SAFD algorithm. It has been demonstrated to 
work in real time on PC type processors and could be implemented 
on such a system at MSFC or other test facilities. 

Validation tests using data from ninety two (92) actual engine tests 
have shown that it is not sensitive to sensor failures and will not 
cause inadvertent shutdowns during nominal engine runs. However, 
because it is a model driven algorithm, it is engine specific; and 
changes to the basic operating characteristics of the engine would 
require changes to the model. For instance, the SSME with the Pratt 
& Whitney ATD pump cannot be monitored with the existing 
algorithm without modifications. Due to this constraint, it is also 
difficult to utilize older engine test data for validation purposes due 
to changes in hardware configurations and sensor measurements. 

Data from past actual SSME failure tests were used to validate the 
sensitivity of the algorithm to engine failures. Seventeen test cases 
were sufficiently similar in operation to the current configuration to 
provide valid comparisons. Twelve of the failures were detected 
significantly earlier than the existing redlines and sufficiently early 
to limit or prevent major damage. Of the five cases in which the 
algorithm detected the failure concurrently with the redlines, three 
were structural failures with no prior warning and the fourth was 
due to a fire with insufficient data available to determine if any 
delay occurred between the onset of the fire and the redline 
detection. The fifth case did exhibit a data spike about 100 seconds 
prior to shutdown which may be detectable by the algorithm if a 
different data filter is developed. It should be noted that these five 
cases did signal shutdown at the same time as the redlines, but are 
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listed as "none" in the table in Appendix A because they were not 
early detections. 

The RSA has reached the maturity level to provide early failure 
detection during ground tests of the SSME. This method is readily 
adaptable to other engine configurations by use of engine specific 
models and could prevent major incidents and concomitant damage 
from test failures. 


2.0 Introduction 

The Rocketdyne Safety Algorithm was originally conceived on the 
LeRC contract Health Management for Rocket Engines . NAS3-25625, 
and is described in the Final Report, NASA CR 185223 (ref 1), and 
other publications concerning the RSA (references 2 through 4). In 
this study contract, the basic form of an algorithm was developed to 
detect failures in the Space Shuttle Main Engine (SSME) ground tests 
earlier than by other available means. The algorithm utilizes a real 
time model of the SSME to compare measured values of key 
parameters with predicted values from the model. The difference 
between the corresponding measured and predicted values for these 
key parameters are normalized, weighted, and summed. This 
weighted sum is then compared to an empirically determined 

threshold for detection of off nominal engine performance. After 
completion of the contract, a Rocketdyne IR&D program was 
conducted to develop a real time model which could be used in 

implementing the algorithm. An influence coefficient engine model 
was developed and a computer test environment was assembled for 
use in validating the algorithm operation. The current contract was 
begun after the IR&D program had successfully concluded. The 
objective of this program is to develop the algorithm to a form which 
can be tested against previous engine test data sets for validation 
and coded for use with actual engine test data at MSFC. The program 

used the model developed on IR&D and the algorithm concept 

previously developed on contract as starting points. Modification of 
the model for a robust system was the initial effort required. 
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3.0 Algorithm Development and Test 

3.1 Algorithm Development 


This Task was begun with an investigation into handling the engine 
startup transient, the development of sensor validation, and 
preliminary ideas for real time automated test case runs. The engine 
model available at the beginning of this Task was developed on IR&D 
and exhibited larger variations during the transient portions of the 
engine test run profile than during the steady state portions, 
particularly during the start transient. Although, the algorithm 
system was not required to operate during the transient portions for 
this Task, it was a desirable feature. In order to more closely predict 
the engine parameters during the startup transient, the Digital 
Transient Model was examined to see if it could lend any processing 
accuracy to the existing safety model algorithm. The digital transient 
model involves predicting engine parameters based upon valve 
position with no inlet conditions using an Adams-Bashforth 
numerical integration routine. For stability this routine requires a 
maximum cycle time of 5 milliseconds over which to iterate. In 
order to determine if this model would be feasible, the integration 
routine was run on eleven computer systems including SUN, PCs and 
Masscomp. The results showed that of the available computers, only 
the Masscomp 68020 -20MHz and an AT 386-33MHz (with math co- 
processor) could run this model if changes were made to reduce 
processing time. Because of this time constraint and the uncertainty 
that this could be implemented successfully in the available 
hardware within the funding constraints of this program, it was 
decided that the Real Time Safety Model (RTSM) developed under 
Rocketdyne IR&D would be used for this contract. 

In an attempt to improve the RTSM, refined power level curve fit 
coefficients for use during start were generated to model engine 
performance during the first five seconds. The start transient 
produced spikes in the weighted sum until a settling of parameters 
occurred. The initial algorithm paused until 5 seconds had passed 
before starting to model the engine, and did not use actual observed 
sensor data until 17 seconds into the test. As a result, weighted sum 
indication of engine performance prior to 17 seconds exhibited 
greater variances than during the later steady state portions. To 
better monitor engine health during the startup transient, the 
algorithm was revised to utilize coefficients derived from previous 
engine tests. Start transient modeling now begins when the first 
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data point is received (0.0 to 0.2 seconds) and the H2 inlet pressure 
is measured. Additionally, the first point of initialization for the 
model was moved from twenty seconds from engine start to five 
seconds from engine start in order to reduce the weighted sum and 
improve the performance during the first twenty seconds. This 
change effectively corrected the start -up transient spiking problem. 


The dependent parameters monitored in the algorithm are listed 
below: 


YNOM(i,01) 

YNOM(i,02) 

YNOM(i,03) 

YNOM(i,04) 

YNOM(i,05) 

YNOM(i,06) 

YNOM(i,07) 

YNOM(i,08) 

YNOM(i,09) 

YNOM(i,10) 

YNOM(i,ll) 

YNOM(i,12) 


= HP FT Discharge Temp A 
= HPFT Discharge Temp B 
= HPOT Discharge Temp A 
= HPOT Discharge Temp B 
= FPOV Position 
= OPOV Position 
= HPFT Speed 
= HPOT Speed 
= OPB Pc 
= FPB Pc 

= MCC Hot Gas INlet Pressure(HG IN) 
= MCC Pc 


Several ideas for measurement validation were investigated for use 
in the algorithm. Reasonableness tests, sensor data weighting, 

fleetwide operating envelopes, and data trending are all possible 
methods for processing data. 

A real time automated test system was configured using a 386 PC to 
take the place of an existing computer because of the speed and disk 
capacity to handle the one hundred nominal SSME tests to be input to 
the model. This PC took the validated, averaged and formatted SSME 
data and processed the model to generate the weighted sum as an 

indication of engine health. 

Execution of the tests of the algorithm was automated to minimize 
time and labor associated with evaluation of the large number of 

hotfire data sets. The existing SSME hot-fire data sets were 

formatted to provide inputs compatible with the Real Time Safety 
Model software. In addition, an automated routine was prepared 
which accessed each hotfire data set in sequence, executed the safety 
algorithm in real time using that data set, and stored a time plot of 
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the algorithm's output (i.e. weighted sum) into an output datafile 
which could be reviewed off-line. 

A 386 PC was integrated and tested in the real time model system to 
provide data storage for both the input SSME hot-fire data files and 
the model's weighted sum output data files. This replaced the slower 
286 machine previously used in the development phase. 

Other work was performed to improve the effectiveness of the safety 
algorithm in preparation for testing. The first initialization point was 
moved to 10 seconds to help insure a smooth transition from engine 
startup onwards. A routine for two point initialization to match 
model parameter curves with measured data at 65% and 100% power 
levels was coded and tested. 

The standard deviations between measured and estimated data for 
the start transient were originally higher than those for mainstage 
data. In order to address this problem, two different arrays for start 
and mainstage calculation of the weighted sum were coded, which 
resulted in significantly better algorithm results when run on non- 
failure test data. This enabled the algorithm to be operational from 
engine start to shutdown. 

Development of the necessary interface for operation on SAFD 

hardware was performed. Parameter calling, timing and receipt of 
parameters were established for the algorithm. Interface 
specifications for communicating with the SAFD system were 
received from Huntsville and reviewed. The specifications were 
reviewed to organize and identify what modifications must be made 
to the model to run on the SAFD platform such as how often the RSA 
algorithm is to be invoked, which CADS and which Facility data are 
needed for the RSA algorithm, etc. The SAFD system already has 
upper and lower limits established for the sensors. Additional limit 
checks may not be needed if they are in agreement with the 

established limits. An algorithm map was created which lists engine 
parameters and information and how often they are accessed. 

A modification was made to the initialization procedure of the 

algorithm which increased its sensitivity. The algorithm now 

initializes each of the twenty two sensors used against the model 
rather than treating the redundant sensors as one and initializing 

only the resultant ten parameters. By this method, each PID is 

initialized to its behavior. This increased the sensitivity because the 
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redundant sensors normally have a small deviation from one 
another. 

A problem caused by an artifact of the test stand calibration 
procedure was investigated and eliminated. After a pre-test sensor 
checkout, one of the channels of the HPFT speed sensor occasionally 
does not clear, leaving a false indication of speed. The code was 
changed so that, at initialization, it checks to see if the HPFT speed is 
less than 4000 RPM for the first 0.7 seconds. If it is, the algorithm 
proceeds normally; otherwise, the speed is not used until real data 
from the engine is received. 

The primary method developed for "soft" sensor failure detection 
involved a "second deviation" of the sensor value. In this process, 
the deviation from the predicted model value of each parameter is 
compared with the corresponding average deviation of all other 
correlated engine parameters to determine if a parameter is showing 
a deviation much different from the others. This "second deviation" 
is used to weight the value of each parameter to the overall sum. 
The weighted sum calculation was reviewed, changed and coded for 
this "soft" sensor failure detection. For each parameter, the deviation 
from the predicted model value is subtracted from the average of all 
other parameter deviations; and the absolute value of the result 
becomes the argument to the hyperbolic secant function which 
returns a value between 0 and 1. This function result is used as a 
non-linear weight for the parameter deviation in the weighted sum 
calculation. The value of a specific parameter deviation will decrease 
in importance as it varies from the average of the other sensors. The 
code includes an adjustable slope for the hyperbolic secant function 
to change performance. In addition, strategies to place sensors into 
related subgroups in which sensors could be compared on a 
correlated basis were developed. 

A Sensor Validation module was developed based on the concept 
that the hot gas flow path on both the fuel and oxidizer systems is 
considered to be a continuous series system beginning at the oxidizer 
preburner valve (OPOV/FPOV) and ending at the main combustion 
chamber pressure (MCC Pc). For example, if the OPOV were to begin 
to close, then the environments at all locations following in the flow 
path would drop. If the OPOV appears to be below the nominal 
position, and the OPB reads at, or above, nominal, then one of these 
two sensors is wrong. Similarly for the case of the OPB reading 
below nominal and the HPOP speed is at or above nominal. 
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The validation begins by first calculating the average of the vnd’s 
(validated normalized deviations). This is done separately for the 
fuel and LOX sides. Each sensor's contribution is removed from the 
average, then its vnd is compared to the average of the remaining 
vnd’s. The algorithm then computes vnd’s for fuel, LOX, and for a 
total. As each sensor is considered for validation, its contribution to 
the vnd is subtracted from the average. Upstream sensor deviations 
that read lower than downstream averages are declared invalid. If 
MCC deviations are higher than the upstream average, then those 
sensors are declared invalid. 

3.2 Validation Testing 

Ninety of 100 nominal tests selected for this contract were run 
against the Rocketdyne Safety Model. Most of the ninety tests in the 
initial runs produced weighted sums below thirty throughout the 
engine test. A few had sums between thirty and fifty. Because this 
larger weighted sum occurs after power level changes, it was 
surmised that the algorithm did not wait long enough before 
returning to using the commanded MCC Pc rather than the measured 
MCC Pc [during power level changes the algorithm operates on the 
measured MCC Pc rather than the commanded MCC Pc]. This higher 
WSUM was not considered a problem since the failures previously 
observed have weighted sum approaching one hundred, but a more 
consistent value for non-failure cases is desirable in order to detect 
failures at the earliest possible time. 

Problems were encountered in attempts to run the algorithm with 
the failure test cases utilized in the SAFD study. The first difficulty 
was due to a problem with the algorithm which was corrected by 
modifying the code. The weighted sum calculation was changed to 
limit the contribution to the weighted sum so that no single 

parameter could cause wsum to indicate failure. This was done by 
setting the maximum npd (npd is the normalized parameter 
deviation) = 4. The weighted sum equation was revised and is now 
equal to l+npd^ for each parameter so that the maximum 

contribution by a single parameter is 65. If no PIDs are available for 

a particular parameter then its contribution to wsum is 1 (which is 

consistent with normalized parameter deviation = 0). Several test 
runs were performed with the modified algorithm to confirm that it 
operated as desired, but continued problems were encountered. The 
wsum general profile for the SAFD failure case data tends to run at a 
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higher level and with spikes and level changes not seen in nominal 
test cases for current tests. This was traced to the fact that the 
engine configuration has changed in significant ways since the SAFD 
failures occurred and the current engine model does not accurately 
predict the dependent parameters in all cases. A similar profile is 
encountered when the algorithm is run with Alternate Turbopump 
engine test data. Failure cases from the current engine configuration 
have the same general value and shape as nominal cases. 
Additionally, many parameters are missing from the old failure cases 
and substitute parameters must be identified and inserted. This has 
not always been successful. As a result, the old failure cases cannot 
be run and evaluated in a straight forward manner. Replacement 
PID's (parameter identification numbers) were selected for 
missing/bad PID's for use in the processing of the engine failure 
tests. The results of several of the runs showed successful early 
failure detections even with the cutoff value set at 60 which is 
required by the higher nominal level of the old cases. 

Due to the algorithm modification, the failure cases and the nominal 
test cases had to be rerun. Several test failure runs were performed 
with the modified algorithm In spite of the difficulties encountered 
with the old test data, many of the failures are indicated by a rise in 
the average value of the algorithm sum. In order to rerun the 
nominal cases, the test data had to be transferred from data tapes to 
the test platform and verified to be in synchronization. Ninety-nine 
nominal tests were run through the Rocketdyne Safety Algorithm, 
but twelve of the nominal tests exhibited spiking in their weighted 
sum after the start and about six of these were of sufficient 
magnitude to cause a test termination. The cause of this problem 
was unknown, and a plan for isolating and correcting the anomaly 
was developed. 

Upon investigation of the anomalous tests, one problem was 
discovered and corrected in the coding which caused one of the 
erroneous test cuts. Upon further investigation, it was determined 
that the remaining "anomalous cuts" were legitimate due to 
anomalous conditions during the tests. These tests were planned for 
"off nominal" conditions or utilized modified hardware and so were 
not terminated early or listed as abnormal. If a "normal" engine 
were to perform in the manner seen in any of these cases during 
ground testing, it would be imperative to shut it down to determine 
what was wrong. 
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After the final correction was made to the coding, the algorithm was 
once more applied to 17 failure cases and the remaining 92 nominal 
cases. Insufficient funding remained to obtain more nominal test 
data cases and prepare them for test in order to perform the 
originally desired 100 nominal test cases. The algorithm worked 
well on the nominal cases and the failure cases. (A summary of the 

nominal and failure cases and the plots of the algorithm value vs. 
time for each test case are presented in Appendix A.) Examination of 
the nominal cases showed that a two level threshold limit could 
effectively be applied. The originally planned limit of 60 worked 
satisfactorily for all cases, but it was apparent that this value could 
be lowered to 40 after the initial 20 seconds of testing. This is not 
unexpected since the model predicts actual performance better after 
the initial calibration checks. This lower limit provides the potential 
for earlier anomaly detection. When applied to the failure cases, the 
algorithm detected the anomalies significantly earlier than the 
failures in 12 of the 17 cases. Of the five tests failures which were 
not detected earlier than the existing redlines, three were structural 
failures with no forewarnings. One was due to a fire in the HPOTP, 
but data are not available on the speed of propagation or what 
triggered the controller to shutdown. The fifth is detectable in a 
spike over 100 seconds early if the current data filtering technique is 
not applied, which indicates that an improvement in the filtering 
method could lead to improved sensitivity. The results of the tests 
indicate that the algorithm is very successful in detecting failures 
significantly earlier than with traditional methods. 

The "SAFD version" of this algorithm was successfully tested on a Sun 
sparcstation at Rocketdyne. A disk with this version and the PC 
version were provided to HSL to verify that the algorithm interfaces 
correctly and operates within expected execution speed for the SAFD 
platform. Tests of the algorithm on the SAFD hardware indicated 
that it works satisfactorily; but due to SAFD system operating time, it 
may not be able to work concurrently with the SAFD algorithm. This 
cannot be verified unless actual tests are performed. It will work 
alone on the SAFD hardware or on a dedicated PC at the facility. A 
future increase in the operating speed of the SAFD hardware would 
also enable it to run concurrently with the SAFD algorithm. The RSA 
code listing is presented in Appendix B. 
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4.0 ALGORITHM DESCRIPTIO N 


The Rocketdyne Safety Algorithm is presented here in a logic 
flowpath. Section 4.1 will describe some of the features in the logic 
diagrams and the general flow of calculations, section 4.2 presents 
the logic diagrams, and section 4.3 contains the values of variables 
used in the algorithm including the influence coefficients. 

4.1 General Flow of Calculations 


• The modules outlined in this document and used in the logic 
flow diagrams only provide logical separation according to 
the algorithm specifications. 

• Using a predetermined array of influence coefficients and 
initial nominal values, estimates of the dependent variables 
are derived from polynomial equations 

• Consistency of redundant measurements is tested against 
maximum and minimum limits 

• Anomalous measurements are deleted from final processing 
or limited in their contribution to the weighted sum 
calculations 


Variables 
ADT[ ] 

COEF[i][m][n] 

D[m] 

DELTAX[n] 
variable values 


Algorithm Data Table value (measured) 

Influence coefficients 
"Normal" parameter deviations 

% deltas between nominal & measured independent 


DELTAY[m] 

DEPNOM[m] 

DMAX 

EST[m] 

INDNOM[n] 

INDVAL[n] 

INFLU [m][n] 

dependent 

ND 

NPD[c] 

PCHK(j] 

P L 


The sum of INFLU(m,n) - DELTAX(n) 

Nominal dependent values 

Maximum number of standard deviations allowed 
between "good" measured and estimated values 
Estimated nominal engine values 
Nominal independent values 
Measured independent parameter values 

Independent parameter influences on each 

parameter 

Normalized measurement deviations 
Normalized parameter deviations 
Maximum delta between the two closest values 
Normalized power level from MCC pressure 


reference 


V[m] 


Array of flags denoting validity of parameters 
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VAL[j] Average of two closest values from ADT for sets of 

3 measurements 

VMAX[b] Values for one standard deviation between deltas 

of normalized, redundant measurements(for 
consistency check) 

VND[] Validated, normalized deviations of each 

measurement (O thru 21 ) 

XNOM Model independent parameter initial conditions 

YNOM Model dependent parameter initial conditions 

Indices 

i = array index (0-3) 

j = subscript for sets of three independent variable measurements 
m = dependent variable subscript (O - 21 ) 
n = independent variable subscript (O - 6) 

a = subscript for redundant pairs of dependent parameters (0-10) 
c = normalized parameter deviation subscript (O - 10) 

Module C 


This module validates the independent parameters: 


LH2 Flow 

LH2 Inlet Pressure 
LH2 Inlet Temperature 
LH2 repressurization flow 


LOX How 

LOX Inlet Pressure 
LOX Inlet Temperature 
LOX repressurization flow 


For these parameters, this module will determine if the sensors for a 
particular measurement agree within defined limits. There are three 
values for each parameter. The two closest values are selected and 
averaged to obtain the validated model input for each parameter. If 
the two closest values differ by more than a predetermined amount, 
that variable is set to zero. This has the effect of disqualifying that 
parameter. The averaged, validated parameters are entered into a 
new array called INDVAL[n]. 


Module D 


By using the independent parameter influences, the 
dependent variables can be estimated more accurately than from 
thrust alone With power level dependence taken into account, 
deviations of a particular measurement from expected levels will 
probably be detected earlier and with greater confidence. Influence 



coefficients can be adjusted as more operational experience is 
accumulated. 


A series of calculations in this module lead to estimated values for 
the 11 redundant pairs of dependent parameters: 

HPFTP Turbine Discharge Temp 

HPOTP Turbine Discharge Temp 

FPOV Position 

OPOV Position 

HPFTP Shaft Speed 

HPOTP Shaft Speed 

OPB Chamber Pressure 

FPB Chamber Pressure 

MCC Hot Gas Injection Pressure 

MCC Chamber Pressure, Channel A 

MCC Chamber Pressure, Channel B 


The DEPNOM[m] array and the INDNOM[n] array are calculated. 

Next, the difference between the independent nominal values, 
INDNOM[n], and the independent measured values, INDVAL[n], is 
divided by INDNOM[n]. This results in the DELTAX[n] array which 
are the percent deltas between nominal and measured independent 
values. 

A group of third order power level equations is then used to 
calculate the INFLU[m][n] matrix which are estimates of the percent 
change in each dependent parameter caused by a 1% perturbation in 
each parameter. 

The total influence on each dependent parameter due to the actual 
perturbations in each independent parameter is calculated by 
summing the products of the influ[m][n] matrix and the deltax[n] 
array. These values are saved as the DELTAY[m] array - the percent 
changes in each dependent parameter. 

The final estimates for dependent parameter values, EST[m], are 
determined by adjusting the nominal dependent values, DEPNOM[m], 
by the percent changes due to independent parameter perturbations, 
DELTA Y[m]. The EST[m] array is calculated by multiplying 
DEPNOM[m] by one plus DELTAY[m]. 
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Module F 


This is the final module for eliminating erroneous dependent 
measurements that are out of their expected range when correlation 
with other measurements is not confirmed. Redundant 
measurements that do not agree within specific limits are also 
eliminated from further use or limited in their contribution to the 
weighted sum calculations. 

The output for this module is the validated, normalized deviation of 
each measurement (all redundant measurements included). There 
are twenty-two (22) measurements. 

An array of twenty -two flags is initialized to keep track of which 
measurements are valid. Each flag, V[m], is set to one or zero. If any 
measurement does not exist (ADT[m] = 0.0) or has been found 
unreasonable in previous modules then the corresponding flag, V[m], 
is set to zero for the remainder of the test. 

Next, the normalized measurement deviations are calculated by 
taking the difference between the measurement (from the Algorithm 
Data Table) and the estimate, EST[m], derived in Module D, and 
dividing by the normal parameter deviation values, D[m]. 

The next step is to test the validity of the normalized deviations just 
calculated. First, an array of maximum differences allowed between 
redundant measurement deviations is initialized with predetermined 
values, VMAX[a]. The minimum average deviation of valid 
measurements that indicates an anomaly is set to the predetermined 
value. 

The absolute value of the difference between each redundant pair of 
normalized parameter deviations, l[ND(al) - ND(a2)]l, is calculated 
and compared to three times the standard deviation for this 
difference, DMAX * VMAX[a]. If this difference is less than three 
standard deviations, the redundant parameters are consistent. 

Module G 

In this module, the validated normalized deviations for each 
of the twenty-two (22) measurements are combined to yield the 
final, normalized parameter deviations, NPD[c]: 



NPD[0] HPFTP Discharge Temperature 

NPD[1 ] HPOTP Discharge Temperature 

NPD[2] FPOV Position 

NPD[3] OPOV Position 

NPD[4] HPFTP Shaft Speed 

NP D[5] HPOTP Shaft Speed 

NPD[6] OPB Chamber Pressure 

NPD[7] FPB Chamber Pressure 

NPD[8] MCC Hot Gas Injection Pressure 

NPD[9] MCC Chamber Pressure, Channel A 

NPD[10] MCC Chamber Pressure, Channel B 


Influence Coefficients 

Influence coefficients are an approximation of output from the Power 
Balance Model (the Power Balance Model is a steady state computer 
model of the SSME capable of modeling effects of variations in 
independent parameters [e.g. engine oxidizer inlet pressure] on 
dependent parameters [e.g. HPOTP turbine discharge temperature]). 
They are based upon the assumption that effects on dependent 
parameters of perturbations in independent parameters are linear 
over ranges of interest. Influence coefficients are used in the 
Rocketdyne Safety Algorithm to estimate values for the 10 
dependent parameters listed above. An estimate for each adjusted, 
nominal dependent parameter value, EST[m], is used as a baseline 
value to compare against the measured value of the dependent 
parameter. 


Influence coefficients consist of three components: 

• third order power level curve fit to the nominal value of 
each dependent parameter: 

DepNom = Aq + Aj*PL + A 2 *PL 2 + A3*PL3 

• third order power level curve fit to the nominal value of 
each independent parameter: 

IndNom = B 0 + Bi*PL + B 2 *PL 2 + B 3 *PL3 

• third order power level curve fits to the percent change in a 

dependent parameter due to a 1% perturbation in an 
independent parameter: 
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Influ = Cq + Ci*PL + C 2 *PL 2 + C3*PL3 


Influence coefficients are used as follows: 

• calculate the percent delta between measured and nominal 

independent parameters: 

DELTAX = ( IndMeasured - IndNom)/IndNom 

• calculate the total effect on the dependent parameter due to 

the variation of each independent parameter from its 

nominal value: 

6 

DELTAY = * DELTAXln] 

n = 0 

• adjust the nominal dependent value for the variations in the 

independent parameters: 

EST = DepNom * (1+DELTAY) 

4.2 Logic Flow Diagrams 

Attached below are the logic diagrams for the Rocketdyne Safety 
Algorithm described above. 




initialize start coefficients 
using 0 to 0.2 s data 
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Safety Monitoring Algorithm^ 
Weighted Sum Calculation J 


Module C 


4 


Module D 


i 


Module F 


validate 6 measured independent parameters by 
comparing 3 redundant values 


f 


yes 




no 


is time < 10 seconds? 


1 


calculate 22 nominal dependent 
values as functions of power level 
(PL) using engine start coefficients 


calculate 22 nominal dependent 
values as functions of PL using 
mainstage coefficients 


calculate 6 nominal independent values as functions of PL 


calculate % deltas between 6 nominal & measured independent values 

I 

baseline the % change in each of 22 dependent 
parameters due to 1% perturbations in each of 6 
independent parameters 


estimate the influence on each dependent parameter due to 
actual % deltas betwen the 6 nominal & measured 
independent values ■ 


check for any non-existent dependent measured values 


i 


I 


is time <10 seconds? 


1 


calculate normalized deltas 
between existing measured & 
predicted dependent values 
using engine start normalizing 


factors 



calculate normalized deltas between 
existing measured & predicted 
dependent values using mainstage 
normalizing factors 



Update the running average of the latest 5 
samples for each normalized delta. 


calculate WSUM by summing all validated 
normalized deltas 
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c 


Safety Monitoring Algorithm 
Sensor Validation 


D 


calculate average of the normalized deviation 

— 8, 




i-1 


* 

calculate the 2nd deviation for each parameter 



Calculate average validated normalized deviations (VND) for LOX and fuel sides 
Remove redundant sensor contributions from average 

Compare VND to average of remaining VNDs 



if VND < (fuel average - standard deviation) 
then Sensor « BAD 
else Sensor = GOOD 


if VND < (LOX average - standard deviation) 
then Sensor - BAD 
else Sensor = GOOD 



Weight valid measurements 
using sech( 52nd ) 


T 

calculate WSUM by summing all validated 
normalized deltas 


18 



Module C #1 

Validate Independent Parameters 


Method: Average of two closest of three values 

• Values set to zero rf two closest differ by 

more than a preselected value 

• Define an array of validated, Model Independent 

Variables: INDVAL(O) thru INDVAL(5) 




Module C #2 

Validate Independent Parameters 

Method: Average of two closest of three values 

• Values set to zero if two closest differ by 

more than a preselected value 

• Define an array of validated, Model Independent 

Variables: INDVAL(O) thru INDVAL(7) 



20 




Module D 



Define 

PL = ADT(34) + 3006 


Calculate nominal dependent values, 
DEPNOM(m) 


22 

values 



No 


Yes 

Shaded blocks done 
in-line as required 


Calculate nominal independent 
values, INDNOM(n) 


Calculate % deltas between 
nominal and measured 
variable values, DELTAX(n) 


Calculate INFLU(m.n) polynomials 
using the influence coefficient 
array (defined in INCLUDE file) 



rr 

lDELTAX(n) - 0.0j 

Y 


Calculate total influence of independent 
variables on dependent variables, 
DELTAY(m) 


6 

values 


22 

values 



Exit 


Calculate estimated dependent 
parameters, EST(m) 


22 

values 


DEPNOM(m) * YNOM(O.m) + YNOM(1 ,m) * PL + YNOM(2.m) * PL A 2 + YNOM(3,m) * PL*3 
INDNOM(n) = XNOM(O.n) + XNOM(1 ,n) * PL + XNOM(2,n) * PL A 2 + XNOM (3,n) * PL A 3 
DELTAX(n) = (INDVAL(n) - INDNOM(n)] + INDNOM(n) 

INFLU(m,n) * COEF(0,m,n) + COEF(l.m.n) * PL + COEF(2.m,n)*PL A 2 + COEF(3,m,n)*PL A 3 
DELTAY(m) = I [INFLU(m,n) * DELTAX(n)] 

M) 

EST(m) = DEPNOM(m) * (1 + DELTAY(m)] 
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MODULE F 

Calculate validated, normalized deviations 


Define arrays D(m) and VMAX(a) 

I 


lsV(a1) = 1 no 



22 



Module G 

Define Normalized Parameter Deviations 


Define the normalized parameter 
deviations in accordance with 
the equation given below 


NPD(a) = (IVND(al)l + |VND(a2)|) / (V(al) + V(a2)) 
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4.3 Constants and Influence Coefficients 


Below are the variable constants, data structures and influence 
coefficients used in the algorithm. 


i = 0 

i = 1 

i = 2 

i = 3 

XNOM(i,0) 6.0 

0.0 

0.0 

0.0 

XNOM(i,l) 30.0 

0.0 

0.0 

0.0 

XNOM(i,2) 100.0 

0.0 

0.0 

0.0 

XNOM(i,3) 37.0 

0.0 

0.0 

0.0 

XNOM(i,4) 164.0 

0.0 

0.0 

0.0 

XNOM(i,5) 5.97545E-02 

0.5719 

6.801E-02 

0.0 

XNOM(i,6) 0.202924 

YNOM START 

1.679 

-2.849E-01 

0.0 

YNOM START(i,0) 

0.000 

0.000 

520.0420 

1777.250 

YNOM START(i,l) 

0.000 

0.000 

520.0420 

1777.250 

YNOM START(i,2) 

0.000 

0.000 

-537.0478 

2449.364 

YNOM START(i,3) 

0.000 

0.000 

-537.0478 

2449.364 

YNOM START(i,4) 

0.000 

0.000 

0.000 

0.000 

YNOM START(i,5) 

0.000 

0.000 

0.000 

0.000 

YNOM START(i,6) 

0.000 

0.000 

0.000 

0.000 

YNOM START(i,7) 

0.000 

0.000 

0.000 

0.000 

YNOM START(i,8) 

63325.11 

0.000 

93179.35 

-124648.3 

YNOM START(i,9) 

63325.11 

0.000 

93179.35 

-124648.3 

YNOM START(i,10) 

0.000 

0.000 

0.000 

0.000 

YNOM START(i,ll) 

0.000 

0.000 

0.000 

0.000 

YNOM START(i,12) 

554.5640 

0.000 

4604.870 

187.0750 

YNOM START(i,13) 

554.5640 

0.000 

4604.870 

187.0750 

YNOM START(i,14) 

1205.108 

0.000 

4952.407 

-1274.236 

YNOM START(i,15) 

1205.108 

0.000 

4952.407 

-1274.236 

YNOM START(i,16) 
0.000 

0.000 

3003.940 

196.7420 

YNOM START(i,17) 

0.000 

0.000 

3003.940 

196.7420 
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YNOM_START(i, 1 8) 0.000 

0.000 

YNOM_START(i, 19) 0.000 

0.000 

YNOM_START(i,20) 0.000 

0.000 

YNOM_START(i,2 1 ) 0.000 

0.000 


YNOM(i,0) 

YNOM(i,l) 

YNOM(i,2) 

YNOM(i,3) 

YNOM(i,4) 

YNOM(i,5) 

YNOM(i,6) 

YNOM(i,7) 

YNOM(i,8) 

YNOM(i,9) 

YNOM(i,10) 

YNOM(i,ll) 

COEF(i,0,l) 

COEF(i,0,2) 

COEF(i,0,3) 

COEF(i,0,4) 

COEF(i,0,5) 

COEF(i,0,6) 

COEF(i,0,7) 

COEF(i,l,l) 

COEF(i,l,2) 

COEF(i,l,3) 

COEF(i,l,4) 

COEF(i,l»5) 

COEF(i,l,6) 

COEF(i,l,7) 

COEF(i,2,l) 

COEF(i,2,2) 

COEF(i,2,3) 

COEF(i,2,4) 

COEF(i,2,5) 

COEF(i,2,6) 

COEF(i,2,7) 

COEF(i,3,l) 

COEF(i,3,2) 

COEF(i,3,3) 

COEF(i,3,4) 

COEF(i,3,5) 

COEF(i,3,6) 

COEF(i,3,7) 


1684.000 

1684.000 

478.000 

478.000 
76.830 
72.590 

17204.000 

4884.000 
-160.000 
-493.000 

132.000 
0.000 

6.7558 10E-02 
-1.70377E-02 
1.339000E-02 
4.452700E-01 
7.812990E-02 
2.482 170E-03 
-6.22494E-04 

6.7558 10E-02 
-1.70377E-02 
1.339000E-02 
4.452700E-01 
7.812990E-02 
2.482170E-03 
-6.22494E-04 

-1.79474E-01 

-9.45672E-04 

-3.15037E-02 

3.09397E-02 

-4.63672E-01 

-5.18391E-03 

1.87290E-03 

-1.79474E-01 

-9.45672E-04 

-3.15037E-02 

3.09397E-02 

-4.63672E-01 

-5.18391E-03 

1.87290E-03 


-746.336 

-746.336 

714.638 

714.638 

-25.4172 

-63.8407 

13012.3 

23191.7 

4604.87 

8396.85 

3003.94 

3006 

3.13551E-01 

1.42833E-02 

-3.66619E-03 

-5.47878E-01 

-3.78262E-01 

-6.37037E-03 

-3.36372E-04 

3.13551E-01 

1.42833E-02 

-3.66619E-03 

-5.47878E-01 

-3.78262E-01 

-6.37037E-03 

-3.36372E-04 

5.309600E+00 

1.754860E-02 

-4.35284E-02 

-2.52552E-01 

1.84179E+00 

1.99238E-02 

3.37052E-03 

5.309600E+00 

1.754860E-02 

-4.35284E-02 

-2.52552E-01 

1.84179E+00 

1.99238E-02 

3.37052E-03 


3006.00 

3006.00 

3006.00 

3006.00 


721.239 

721.239 

110.014 

110.014 

29.1003 

57.6847 

4720.79 

-836.435 

187.075 

-4697.71 

196.742 

0.000 

-6.33715E-01 

-4.504 18E-03 

-3.06278E-03 

3.84827E-01 

1.77979E-01 

2.89615E-03 

3.61408E-04 

-6.337 15E-01 

-4.504 18E-03 

-3.06278E-03 

3.84827E-01 

1.77979E-01 

2.89615E-03 

3.61408E-04 

-1.78281E+00 

-1.05341E-02 

4.25821E-02 

0.00000E+00 

-7.88768E-01 

-8.76068E-03 

-2.36554E-03 

-1.78281E+00 

-1.05341E-02 

4.2582 IE-02 

0.00000E+00 

-7.88768E-01 

-8.76068E-03 

-2.36554E-03 


0.0000 

0.000 

0.000 

0.000 


0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

554.564 

2445.17 

0.000 

0.000 

O.OOOOOE+OO 

0.00000E+00 

0.00000E+00 

O.OOOOOE+OO 

O.OOOOOE+OO 

O.OOOOOE+OO 

O.OOOOOE+OO 

O.OOOOOE+OO 

O.OOOOOE+OO 

O.OOOOOE+OO 

O.OOOOOE+OO 

O.OOOOOE+OO 

O.OOOOOE+OO 

O.OOOOOE+OO 

O.OOOOOE+OO 

O.OOOOOE+OO 

O.OOOOOE+OO 

O.OOOOOE+OO 

O.OOOOOE+OO 

O.OOOOOE+OO 

O.OOOOOE+OO 

O.OOOOOE+OO 

O.OOOOOE+OO 

O.OOOOOE+OO 

O.OOOOOE+OO 

O.OOOOOE+OO 

O.OOOOOE+OO 

O.OOOOOE+OO 
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COEF(i,4,l) 

COEF(i,4,2) 

COEF(i,4,3) 

COEF(i,4,4) 

COEF(i,4,5) 

COEF(i,4,6) 

COEF(i,4,7) 

C0EF(i,5,l) 

COEF(i,5,2) 

COEF(i,5,3) 

COEF(i,5,4) 

COEF(i,5,5) 

COEF(i,5,6) 

COEF(i,5,7) 

COEF(i,6,l) 

COEF(i,6,2) 

COEF(i,6,3) 

COEF(i,6,4) 

COEF(i,6,5) 

COEF(i,6,6) 

COEF(i,6,7) 

COEF(i,7,l) 

COEF(i,7,2) 

COEF(i,7,3) 

COEF(i,7,4) 

COEF(i,7,5) 

COEF(i,7,6) 

COEF(i,7,7) 

COEF(i,8,l) 

COEF(i,8,2) 

COEF(i,8,3) 

COEF(i,8,4) 

COEF(i,8,5) 

COEF(i,8,6) 

COEF(i,8,7) 

COEF(i,9,l) 

COEF(i,9,2) 

COEF(i,9,3) 

COEF(i,9,4) 

COEF(i,9,5) 

COEF(i,9,6) 

COEF(i,9,7) 

COEF(i,10,l) 

COEF(i,10,2) 

COEF(i,10,3) 

COEF(i,10,4) 

COEF(i,10,5) 


-3.27388E+00 

-2.72697E-02 

4.14552E-02 

1.65703E+00 

-7.63844E-02 

-1.04716E-02 

-2.40354E-03 

-3.27388E+00 

-2.72697E-02 

4.14552E-02 

1.65703E+00 

-7.63844E-02 

-1.047 16E-02 

-2.40354E-03 

1.33173E+00 

6.04503E-04 

1.30019E-02 

2.20889E-02 

8.40340E-01 

- 1.191 88E-03 

-9.63561E-04 

1.33173E+00 

6.04503E-04 

1.30019E-02 

2.20889E-02 

8.40340E-01 

-1.19188E-03 

-9.6356 IE-04 

-6.31739E-01 

-1.17272E-02 

1.19841E-03 

3.06835E-01 

-2.97835E-02 

-1.32908E-04 

-1.11418E-04 

-6.31739E-01 

-1.17272E-02 

1.19841E-03 

3.06835E-01 

-2.97835E-02 

-1.32908E-04 

-1.11418E-04 

-1.53861E-01 

7.03175E-05 

-5.88139E-02 

-4.07650E-03 

1.13350E-01 


8.50257E+00 

6.61322E-02 

-8.71833E-02 

-4.27165E+00 

2.17812E-01 

2.78641E-02 

6.10207E-03 

8.50257E+00 

6.61322E-02 

-8.71833E-02 

-4.27165E+00 

2.17812E-01 

2.78641E-02 

6.10207E-03 

-3.68307E+00 

0.00000E+00 

-2.74375E-02 

-6.10283E-02 

-2.36836E+00 

2.29577E-03 

2. 1 3802E-03 

-3.68307E+00 

0.00000E+00 

-2.74375E-02 

-6.10283E-02 

-2.36836E+00 

2.29577E-03 

2.13802E-03 

1.12647E+00 

1.26742E-02 

-4.42708E-03 

-3.73406E-01 

8.02084E-02 

O.OOOOOE+OO 

3.62329E-04 

1.12647E+00 

1.26742E-02 

-4.42708E-03 

-3.73406E-01 

8.02084E-02 

O.OOOOOE+OO 

3.62329E-04 

7.19883E-01 
O.OOOOOE+OO 
9.10559E-02 
7.660 15E-03 
3.59439E-01 


-6.35848E+00 

-4.57526E-02 

5.334435E-02 

2.88507E+00 

-1.79975E-01 

-2.0263 IE-02 

-4.52982E-03 

-6.35848E+00 

-4.57526E-02 

5.334435E-02 

2.88507E+00 

-1.79975E-01 

-2.0263 IE-02 

-4.52982E-03 

3.04542E+00 

O.OOOOOE+OO 

O.OOOOOE+OO 

O.OOOOOE+OO 

1.88330E+00 

O.OOOOOE+OO 

O.OOOOOE+OO 

3.04542E+00 

O.OOOOOE+OO 

O.OOOOOE+OO 

O.OOOOOE+OO 

1.88330E+00 

O.OOOOOE+OO 

O.OOOOOE+OO 

-8.22106E-01 

-5.26788E-03 

2.67544E-03 

2.76127E-01 

-4.00085E-02 

O.OOOOOE+OO 

-2.01927E-04 

-8.22106E-01 

-5.26788E-03 

2.67544E-03 

2.76127E-01 

-4.00085E-02 

O.OOOOOE+OO 

-2.01927E-04 

-3.08034E-01 

O.OOOOOE+OO 

-4.02095E-02 

-5.73094E-03 

-1.50146E-01 


O.OOOOOE+OO 

O.OOOOOE+OO 

O.OOOOOE+OO 

O.OOOOOE+OO 

O.OOOOOE+OO 

O.OOOOOE+OO 

O.OOOOOE+OO 

O.OOOOOE+OO 

O.OOOOOE+OO 

O.OOOOOE+OO 

O.OOOOOE+OO 

O.OOOOOE+OO 

O.OOOOOE+OO 

O.OOOOOE+OO 

O.OOOOOE+OO 

O.OOOOOE+OO 

O.OOOOOE+OO 

O.OOOOOE+OO 

O.OOOOOE+OO 

O.OOOOOE+OO 

O.OOOOOE+OO 

O.OOOOOE+OO 

O.OOOOOE+OO 

O.OOOOOE+OO 

O.OOOOOE+OO 

O.OOOOOE+OO 

O.OOOOOE+OO 

O.OOOOOE+OO 

O.OOOOOE+OO 

O.OOOOOE+OO 

O.OOOOOE+OO 

O.OOOOOE+OO 

O.OOOOOE+OO 

O.OOOOOE+OO 

O.OOOOOE+OO 

O.OOOOOE+OO 

O.OOOOOE+OO 

O.OOOOOE+OO 

O.OOOOOE+OO 

O.OOOOOE+OO 

O.OOOOOE+OO 

O.OOOOOE+OO 

O.OOOOOE+OO 

O.OOOOOE+OO 

O.OOOOOE+OO 

O.OOOOOE+OO 

O.OOOOOE+OO 
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COEF(i,10,6) 

COEF(i,10,7) 

COEF(i,10,l) 

COEF(i,10,2) 

COEF(i,10,3) 

COEF(i,10,4) 

COEF(i,10,5) 

COEF(UO,6) 

COEF(i,10,7) 

COEF(UU) 
COEF(i,l 1 ,2) 
COEF(i,ll,3) 
C0EF(i,ll,4) 
COEF(i,ll,5) 
COEF(i,ll,6) 
COEF(i,ll,7) 

COEF(i,12,l) 

COEF(i,12,2) 

COEF(i,12,3) 

COEF(i,12,4) 

COEF(i,12,5) 

COEF(i,12,6) 

COEF(i,12,7) 

COEF(i,13,l) 

COEF(i,13,2) 

COEF(i,13,3) 

COEF(i,13,4) 

COEF(i,13,5) 

COEF(i,13,6) 

COEF(i,13,7) 

COEF(i,14,l) 

COEF(i,14,2) 

COEF(i,14,3) 

COEF(i,14,4) 

COEF(i,14,5) 

COEF(i,14,6) 

COEF(i,14,7) 

COEF(i,15,l) 

COEF(i,15,2) 

COEF(i,15,3) 

COEF(i,15,4) 

COEF(i,15,5) 

COEF(i,15,6) 

COEF(i,15,7) 

COEF(i,16,l) 

COEF(i,16,2) 

COEF(i,16,3) 


-9.4548 IE-04 
8.458920E-04 

-1.53861E-01 

7.03175E-05 

-5.88139E-02 

-4.07650E-03 

1.13350E-01 

-9.4548 IE-04 

8.458920E-04 

-1.06276E-02 

-3.09662E-04 

-1.11283-02 

7.83933E-03 

-2.95686E-02 

-4.82056E-04 

6.0501 IE-04 

-1.06276E-02 

-3.09662E-04 

-1.11283-02 

7.83933E-03 

-2.95686E-02 

-4.82056E-04 

6.0501 IE-04 

-4.17196E-01 

-3.58206E-03 

1.29300E-03 

1.15425E-01 

-2.13535E-02 

-1.19488E-04 

-8.58883E-05 

-4.17196E-01 

-3.58206E-03 

1.29300E-03 

1.15425E-01 

-2.13535E-02 

-1.19488E-04 

-8.58883E-05 

-2.6668 IE-02 

-2.37199E-04 

-3.22731E-04 

1.00202E-02 

1.17278E-03 

2.65487E-06 

2.85558E-05 

-2.6668 IE-02 
-2.37199E-04 
-3.22731E-04 


4.14947E+00 

O.OOOOOE+OO 

7.19883E-01 

O.OOOOOE+OO 

9.10559E-02 

7.66015E-03 

3.59439E-01 

4.14947E+00 

O.OOOOOE+OO 

4.20930E-02 

-3.51086E-04 

5.70060E-03 

-2.56483E-04 

2.10816-E01 

1.94540E-03 

-1.20095E-04 

4.20930E-02 

-3.51086E-04 

5.70060E-03 

-2.56483E-04 

2.10816-E01 

1.94540E-03 

-1.20095E-04 

3.95345E-01 

2.1991 IE-03 

-4.85527E-03 

-1.52808E-01 

5.46972E-02 

O.OOOOOE+OO 

3.00620E-04 

3.95345E-01 

2.19911E-03 

-4.85527E-03 

-1.52808E-01 

5.46972E-02 

O.OOOOOE+OO 

3.00620E-04 

9.92401E-03 

O.OOOOOE+OO 

O.OOOOOE+OO 

-1.39743E-02 

-3.27962E-05 

-1.72772E-04 

-1.72770E-05 

9.92401E-03 

O.OOOOOE+OO 

O.OOOOOE+OO 


-1.81101E-03 

O.OOOOOE+OO 

-3.08034E-01 

O.OOOOOE+OO 

-4.02095E-02 

-5.73094E-03 

-1.50146E-01 

-1.81101E-03 

O.OOOOOE+OO 

O.OOOOOE+OO 

2.46367E-04 

O.OOOOOE+OO 

1.021 11E-02 

-7.71641E-02 

-7.86290E-04 

O.OOOOOE+OO 

O.OOOOOE+OO 

2.46367E-04 

O.OOOOOE+OO 

1.021 11E-02 

-7.71641E-02 

-7.86290E-04 

O.OOOOOE+OO 

-3.26333E-01 

-6.034 13E-04 

2.59230E-03 

1.16439E-01 

-1.49099E-02 

O.OOOOOE+OO 
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COEF(i,16,4) 
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5.0 Development Elaji 


5.1 Study Plan 

There are several follow on development tasks which should be 
undertaken. These tasks range from further development of the 
algorithm itself to a more complementary integration with the SAFD 
platform. These tasks include, but are not limited to, failure 

detection techniques, improve startup modeling, investigation for 
new dependent parameters, and sensor validation. 


Recommendations 

Despite the success of the Rocketdyne Safety Algorithm, several 
areas of the algorithm need to be addressed and refined. 


•sensor validation 

•accounting for failure of PID 129 or 130 (MCC PC 

A1 & A2) and calculating power level 

•making wsum a function 

•dealing with measured pc during transients 

There are indications that the accuracy of predictions declines as 
engine power level varies from the initialization point. A potential 
solution is to make wsum a function of how far away the current 
power level is from what the power level was during the 
initialization point. The wsum would be scaled down in cases of 
extreme power level variation to account for less accurate 
predictions. 

PIDs 129 and 130 (MCC PC A1 & A2) are used to calculate power 
level during power level transients. There is currently no way to 
replace one or both of these PIDs in the power level calculation if 
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these PIDs were to fail. A methodology to do this should be 
incorporated into the algorithm. 

New Dependent Pa rameters 

Currently 11 pairs of redundant parameter measurements are used 
by the algorithm to determine dependent parameter values. It will 
be a worthwhile task to investigate additional parameters which may 
be useful as dependent parameters. Prime candidates would be 
other PIDs which are redundant measurements since they can be 
used to help validate each other. PIDs which fall into this category 
include: 

LPFTP and LPOTP speeds 

LPFP, HPFP, LPOP, HPOP and PBP discharge pressures 

HPFP and HPOP inlet pressures 

HPFP inlet temperature 

PBP discharge temperature 

MCC oxidizer injection temperature 

New dependent parameters should not be limited to redundant PIDs, 
however. Several single measurement PIDs are very indicative of 
engine condition. The addition of these PIDs to the set of dependent 
parameters should be investigated as well. 

Sensor Validation 

For a full treatment of sensor validation, the algorithm should not 
only look at multiple sensors and the comparison of their individual 
values, but also the interaction of different sensors. This sensor 
validation would involve looking at flow paths through the engine 
and determining the series of sensors which are linked in behavior 
so that it can be determined if an individual sensor has gone awry. 
For example, the hot gas flow path on both the fuel and oxidizer 
systems is considered to be a continuous series system beginning at 
the oxidizer preburner valve (OPOV/FPOV) and ending at the main 
combustion chamber pressure (MCC Pc). If the OPOV were to be 
closed off, then the environments at all locations following in the 
flow path would drop to zero. It follows that if, instead of closing it 
completely, the OPOV position were but reduced, then the following 
environments would be correspondingly reduced. If the OPOV 
appears to be below the nominal position, and the OPB reads at or 
above nominal, then one of these two sensors is wrong. Similarly for 
the case of the OPB reading below nominal and the HPOP speed at or 
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above nominal. This is an example of a rule that would be developed 
and coded into the algorithm. 

Startup modeling 

The current power level based curves used to predict parameter 
values during engine startup are fits to data from a very limited 
number of tests. Three areas have been identified as having 
potential for improvement in startup modeling. 

1) Use the DTM to generate nominal SSME startup data, and generate 
new curve fits from these data. The resulting curves should be more 
representative of a nominal engine. 

2) Base at least part of the startup modeling on commanded positions 
of propellant valves. This should improve the prediction of MCC and 
preburner primes and other start fluctuations. 

3) Use start envelopes to define upper and lower bounds of 
parameter values. The envelopes will be wider at certain operating 
times when confidence in predictions declines. 

Engine Diagnostics 

Three areas where the diagnostic capabilities of the algorithm could 
be further developed are the baseline WSUM, data trending and 
failure pattern recognition. 

In the current code, an engine shutdown command is issued when 
the wsum value exceeds a threshold value. This value could be the 
result of numerous parameters deviating from nominal operation or 
just a few parameters. It would be beneficial if there were a "yellow 
line" or region wherein the weighted sum value was above that of 
nominal engine operation but below the shutdown limit. In this 
manner, the engineers would know that the engine is performing out 
of norm but not in a manner representing critical failure. 

Failure detection exists in the Rocketdyne Safety Algorithm as a flag 
indicating failure when the weighted sum value exceeds the 
threshold value. This is the simplest way to detect failure and in the 
case of catastrophic failure should command engine shutdown before 
engine damage. However, a more sophisticated approach would 
involve looking at the data for trends. For example, a slope-average 
algorithm might be employed to detect slow failures. In this case, 
the slope of the averaged data would be examined rather than the 
averaged data. In this manner, the RSA would be more sensitive to 
subtle changes in slope and could track increasing slope-averages for 
trends which could indicate engine failure. 
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The failure pattern recognition task would involve studying the 
details of all the dependent parameters and their variations for 
various failures. 

Algorithm Validation ... . 

Once new dependent parameters, sensor validation, improved start- 
up modeling and engine diagnostics have been incorporated into the 
this algorithm, it must be validated by monitoring nominal tests and 
a variety of types of failure tests. This task would involve running 
the algorithm with previously transferred tests and transferring the 
more recent SSME tests available to Rocketdyne. 

Feasibility Study of Pratt & Whitnev Pump 

This sub task would involve studying and evaluating the P&W pump 
performance data to determine a methodology for incorporating the 
P&W pump into the model. 

Work Breakdown Structure 

The proposed effort is broken down into the following WBS elements: 

1.0 improve Algorithm Robustness 

1.1 Develop new dependent parameters 

1.2 Develop more robust sensor validation methodology 

1.3 Improve engine startup model 

2.0 Develop engine diagnostics module 

2.1 Determine baseline WSUM 

2.2 Develop data trending methodology 

2.3 Develop failure pattern recognition module 

3.0 Perform P&W pump feasibility study 

4.0 Incorporate changes into RSA 

4.1 Modify algorithm interface & code 

4.2 Perform algorithm validation testing 

5.0 Task Management 

5.1 Program management 

5.2 Travel 

5.3 Reporting 

The cost of the complete program described above including travel, 
manpower, and material is estimated at $250,000. The logic flow 
diagrams for the specific tasks and a schedule follow. 
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IMPROVE STARTUP MODELLING 
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APPENDIX A 

RSA VALIDATION TEST RESULTS 
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FAILURE TEST LIST 


WSUM threshold 60 from 0-20 Seconds, 40 from 20 Seconds - C/O 


TEST# 

CUT-OFF 

[sec] 

RSA 

CUT-OFF 

COMMENTS 

A1331 

233 

NONE 

Injector Lox post failure - very rapid - no time to respond prior 
to C/O 

A1307 

75 

22 

RSA C/O at 22 may be due to model differences causing level 
shift - if so RSA /C/O occurs between 50 and 64 seconds 
depending upon amount of shift 

A3259 

101 

20 

RSA C/0@ 20 to 22 seconds even if model difference shift 
has occurred. 

A3 175 

115 

NONE 

Duct structural failure - no early warning for any system to 
detect using existing sensors 

A1225 

255.63 

145 

Lox ignition due to valve fretting resulting from high vibrations - 
RSA detected cause which eventually led to ignition 

| 

A1110 

74 

NONE ; 

Fire occurred in HPOTP - C/O by controller - don’t have data 
on speed of failure propagation or redline trigger 

A1340 

405 

NONE 

Indicators of failure occur at 290 seconds, but current 
algorithm does not detect due to sharpness of data spike - 
without filter, the failure is detected. Future work may provide 
improved filter technique which would detect this failure 

A1363 

250 

135 

j 

HPFTP Turnaround duct cracked - RSA algorithm may have 
detected as early as 30 seconds, but certainly at 1 35 seconds 

A1436 

611 

NONE 

Coolant liner buckled at 61 0.36 - no time to react 

A1364 

392 

40(25) 

Kaiser hat nut failure early in test propogates to worse and 
worse failures as test continues until engine is destroyed. First 
RSA failure signal is about 40 seconds if model shift is 
assumed (25 seconds, if not) - Algorithm shows continued 
buildup of distress until cutoff occurs 
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FAILURE TEST LIST 


WSUM threshold 60 from 0-20 Seconds, 40 from 20 Seconds - C/O 


TEST# 

CUT-OFF 

[sec] 

RSA 

CUT-OFF 

COMMENTS 

A2209 

823 

595 

HPFT end dome lock tab breaks loose and propagates 
damage - not catastrophic 

A2249 

450 

150 

Initial damage began early in test & slowly propogated. HPFP 
pump cavitation began at 108 seconds and slowly increased. 
KelF rub ring fails at 374 seconds, first stage turbine blade fails 
at 450 seconds. RSA detected failure during early cavitation 
period - algorithm shows increase in distress throughout test 

A2095 

51 

8 

Shutdown due to preburner pump radial accel redline - eroded 
turbine damage and 8 main injector Lox posts eroded • RSA 
detects early indication at 8 seconds, and shows buildup of 
distress beginning at 30 seconds 

A1346 

500 

20 

Damage to HPFT and MCC liner - Strong indication of 
problem by RSA at about 20 seconds 

A1362 

500 

240(20) 

HPOT and HPFT - damaged - RSA nominal value appears to 
be shifted up due to model difference - failure clearly indicated 
at 240 seconds even if shifted - if not shifted, failure noted at 
20 seconds 

All 83 

50 

2 

Injector face burn through due to structural failure - self limiting 
- RSA indicates major problem between 2 and 10 seconds 

A2428 

204.1 

130 

Good example of recent failure (not SAFD database ) - RSA 
shows clearly distress buildup as time progresses with cutoff 
at 1 30 seconds 
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APPENDIX B 
RSA CODE 




# include <math.h> 


/************************** 


ROCKETDYNE SAFETY ALGORITHM 
ALGORITHM MOOULE INITO 


PURPOSE: 


This module initializes the local variables of the RSA module. It is 
called prior to processing a test with RSA. The SAFD platform requires 
a separate initO routine in addition to the RSA main module, thus 
this code. 


This module does nothing except tell RSA to initialize itself. 


r****************** *********************************************************' 


/ 


*********************** 


ROCKETDYNE SAFETY ALGORITHM 
ALGORITHM MOOULE RSA() 


PURPOSE: 


This module contains the code for the RSA. The interface to this module 
from the outside world is the parameter list ONLY. 



/** RSA FLOW-CONTROL VARIABLES **/ 
/** (not directly used for calculation of redlining sun) **/ 
/** (Used to track power levels, keep track of time, etc.) **/ 


/** RSA FLOW-CONTROL VARIABLES **/ 
/** (not directly used for calculation of redlining sun) **/ 
/** (Used to track power levels, keep track of time, etc.) **/ 
/** **/ 


/* Nuneric constants needed for data structure sizes: */ 


#def ine MAX 80 
#def ine ARRAYMAX 2500 
#def ine maxadt 100 
#def ine maxest 22 
#def ine maxdepnom 22 


/♦ constants for data structure sizes */ 
/* Mostly used in USUMO routine. V 
/* and support sdbs (engine model) */ 


/* Sampling is done to average data for the power level signature. V 
/* This signature is otherwise known as the “initialization point" */ 
/* These variables define sampling times in seconds: V 


PKfcCSDlttS PAGE BLANK MOT FILMED 
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# include <math.h> 


ROCICETDYNE SAFETY ALGORITHM 
ALGORITHM MOOULE 1NITO 

PURPOSE: 

This module initializes the local variables of the RSA module. It is 
called prior to processing a test with RSA. The SAFD platform requires 
a separate init() routine in addition to the RSA main module, thus 
this code. 

This module does nothing except tell RSA to initialize itself. 



ROCICETOYME SAFETY ALGORITHM 
ALGORITHM MOOULE R$A<) 

PURPOSE: 

This module contains the code for the RSA. The interface to this module 
from the outside world is the parameter list ONLY. 



/* Numeric constants needed for data structure sizes: */ 


♦define MAX 80 /* constants for data structure sizes */ 

♦define ARRAYJNAX 2500 /+ Mostly used in USUMO routine. V 

♦define max_adt 100 /* and support subs (engine model) */ 

♦define max_est 22 
♦define maxjdepnoei 22 

/* Sampling is done to average data for the power level signature. V 
/* This signature is otherwise known as the "initialization point" V 
/* These variables define sampling times in seconds: V 
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static float time_that_we_begin_sampl ing = 5.9; 
static float time_main__stage_is_r cached = 5.9; 
static float t ime_to_wai t = 4.0; 

static float endofsampl ing_period = -1; /* set at run-time */ 

/* variables to make us wait after power level transients so that we V 
/* don't model incorrectly. 01/16/92 — D'Valentine/RVdH V 

static float transient_count = 0; 
static float transient^ lag = 1; 

/* Arrays used for sampling, averaging of data to produce engine */ 

/* signatures: V 

static float VH[22], woax_tot, wfax_tot, deltayj. 1223 ; 

static float averageji [max_adt] ; /* averages of data for power level */ 

static float averagej. [max_adt] ; /* modelling. */ 

/* Counters and flags to control sampling process: V 

static int sample_count = 0; /* nunber of samples taken so far V 

static float Pc_old = -999999.0; /* previous history memory NccPc V 

static int donees amp l ing = 0; /* Are we done sampling? */ 

static int very_f irst_data _point = 1; /* for initial processing... */ 



/* INPUT/OUTPUT PARAMETER LIST TYPES for calling the RSA module, */ 
/* done the stupidest (but most transportable) way possible, V 

/* as the specifications require us to do V 

/* This code is identical in both RSA.C and SAFD.C, and really V 
/* should be extracted to a separate include file for generality V 


struct inp { float adtO, adtl, adt2, adt3, adt4, adtS, 

adt6, adt7, adt8, adt9, adtIO, adtll, 
adt12, adt13, adt14, adt15, adt16, adt17, 
adt18, adt19, adt20, adt21, adt22, adt34, 
adt38, adt39, adt40, adt41 , adt42, adt43, 
adt47, adt48 f adt49, adtSO, adt51 # adt52, 
adt53 f adt54, adt55, adt56, adt57, adt90; 


long int enginestatus; >; 


struct inp tags { long int adtO, adtl, adt2, adt3, adt4, adt5, 

adt6, adt7, adt8, adt9, adtIO, adtll, 
adt12, adt13, odt14, adtIS, adt16, adt17, 
adt18, adt19, adt20, adt21, adt22, adt34, 
adt38, adt39, adt40, adt41, adt42, adt43, 
adt47, adt46, adt49, adtSO, adt51. adt52, 
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adt53, adtS4, adt55, adt56, adt57, adt90 f 
engine_status; 

>; 


struct outp i long int shutdown; 

float adtO, adtl, adt2, adt3 f adt4, adt5, 

adt6, adt7, adt8, adt9 r adtIO, adtll, 
adt12, adt13, adt14, adt15, adt16, »dt17 # 
adt18, adt19, adt20, adt21, adt22, adt34, 
adt38, adt39, adt40, adt41 # adt42, adt43, 
adt47, adt48, adt49, adt50, adt51, adt52, 
adt53, adt54, adt55, adt56, adt57, adt90, 
estO, estl, est2, est3, est4, est5, est6, 
est7, est8, est9, estIO, estll, est12 r 
est13, est14, est15, est16, est17, est18, 
est19 # est20, est21, wsum; 

>; 


struct outp_tags C long int shutdown, adtO, adtl, adt2, adt3, adt4, adt5, 

adt6, adt7, adt8 # adt9, adtIO, adtll, 
adtl 2, adt13, adt14, adtIS, adt16, adt17, 
adt18, adt19, adt20, adt21, adt22, adt34, 
adt38, adt39, adt40, adt41, adt42, adt43, 
adt47, adt48, adt49, adtSO, adtSI, adt52, 
adtS3, adt54, adt55, adt56, adt57, adt90, 
estO, estl, est2, est3, est4, est5, est6, 
est7, est8, est9, cstIO, estll, cst12, 
est13, est14, est15, est16, est17, est18, 
est19, est20, est21, wsun; 

>; 


/* ACTUAL POINTERS TO INPUT/OUTPUT ARRAYS FOR GLOBAL USAGE: */ 

struct inp *indata; 
struct outp *outdata; 
struct outp_tags *out_tags; 

/* Local Oata Table to which data is copied from outside sources: V 
static float ADTC100]; 


static float valid]; /* misc. variables used for numeric calc. */ 

static float nval; /* Mostly used in either USUMC) or wfax() */ 

static double indnom[7]; /* or woax() routines. V 

static float w[123 ; 

static double depnom[22]; 

static double dep_star[22] ; 

static double influ [22H7] ; 

static float indval[7]; 

static float vnd[22]; 

static float v[22]; 

static float npd[12]; 

static float ndt22]; 

static float ed_nd[22] ; 

static float f ilt_nd[22] [10]; 

static float ad; 
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static float deltay[22]; 
static float del tax (7]; 

static float PCHKH8) = <180.0, 10.0, 25.0, 6.0, 18.0, 85.0, 1.0, 1.0); 


/* The Matrix representing our models of engine parameters (derived later): */ 

static float est[22] = < 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 

0 . 0 , 0 . 0 , 0 . 0 , 0 . 0 , 0 . 0 , 0 . 0 , 

0 . 0 , 0 . 0 , 0 . 0 , 0 . 0 , 0 . 0 , 0 . 0 , 

0 . 0 , 0 . 0 , 0 . 0 , 0.0 ); 

/* matrixes added to enhance ability to handle power level changes */ 

/* Added 9/19/91 by Mark O'Valentine. */ 

static float d_st[22] = < 141.48, 134.71, 167.83, 160.74, *1.000, *1.000, 

-1.000, -1.000, 5252.1, 5251.9, -1.000, -1.000, 

688.65, 688.65, 176.85, 161.71, 80.222, 48.221, 

5.7192, 5.9477, 25.349, 26.213 ); 

static float djns [22] = < 63.614, 66.340, 112.44, 104.73, 4.0549, 4.0648, 

2.6113, 2.5301, 302.95, 306.82, 424.91, 424.91, 

95.727, 95.727, 194.03, 191.24, 30.237, 22.150, 

10.200, 10.790, 11.009, 12.525 ); 

static float vmax_st[11] - < 0.2144, 0.1708, 0.0000, 0.0000, 0.0150, 0.0000, 

0.0000, 0.9965, 2.5076, 2.0956, 0.5132 >; 

static float vmax_ms[11] = < 1.7390, 0.7460, 0.1050, 0.2639, 0.5132, 0.5000, 

0.5000, 0.3350, 2.4651, 0.7917, 0.8479 ); 


/* Cross-reference arrays used for automating a clunsy numbering scheme */ 

/* that we use to identify ADT/EST parameters, thus the repetitions of */ 

/* 7,7,8,8,9,9... etc. RVdH 2/3/92 */ 

static short vmaxcross til] = <0,1 ,2, 3, 4, 5, 6, 7, 8, 9, 9); 
static short vndcross [11] » <0,2,4,6,8,10,12,14,16,18,20); 

static short ndcross [22] = <0,1,2,3,4,4,5,5,6,6,7,7,8,8,9,9,10,10,11,11,11,11); 
static short crossref [6] [3] = <35, 36, 37, 

38, 39. 40, 

47, 48, 49. 

41, 42. 43, 

50. 51. 52, 

44, 45, 46 >; 

/* Matrixes of coefficients to use in polynomial engine model. Used in */ 

/* WSUMO and support routines which set up model. */ 

static float XNOM [71 [4] = < 6.0, 0.0, 0.0, 0.0, 

30.0, 0.0, 0.0, 0.0, 

100 . 0 , 0 . 0 , 0 . 0 , 0 . 0 , 

37.0, 0.0, 0.0, 0.0, 

164.0, 0.0, 0.0, 0.0, 

5.97545E-02, 0.5719, 6.801E-02, 0.0, 

0.202924, 1.679, -0.2849, 0.0); 

static float ED_N0H[7]; 

static double ynomt[22]; /* need not be initialized; derived. V 

static float Y1NI [22] [4] ; 

static float YN0M [22] [4] = < 1684.00, -746.336, 721.239, 0.0, 

1684.00, -746.336, 721.239, 0.0, 
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478.000, 

478.000, 
76.83, 
76.83, 
72.59, 
72.59, 

17204.0, 

17204.0, 

4884.0, 

4884.0, 
-160.00, 
-160.00, 
-493.00, 
-493.00, 

132.00, 

132.00, 

0 . 0 , 

0 . 0 , 

o.o, 

0 . 0 , 


714.638, 

714.638, 

-25.4172, 

‘25.4172, 

-63.8407, 

-63.8407, 

13012.3, 

13012.3, 

23191.3, 

23191.3, 
4604.87, 
4604.87, 
8396.85, 
8396.85, 
3003.94, 
3003.94, 

3006.0, 

3006.0, 

3006.0, 

3006.0, 


110.014, 

110.014, 

29.1003, 

29.1003, 

57.6847, 

57.6847, 

4720.79, 

4720.79, 

-836.435, 

-836.435, 

187.075, 

187.075, 

-4697.71, 

-4697.71, 

196.742, 

196.742, 

0 . 0 , 

0 . 0 , 

o.o, 

0 . 0 , 


0 . 0 , 

0 . 0 , 

0 . 0 , 

0 . 0 , 

0 . 0 , 

0 . 0 , 

0 . 0 , 

0 . 0 , 

0 . 0 , 

0.0, 

554.564, 

554.564, 

2445.17, 

2445.17, 

0 . 0 , 

0 . 0 , 

0 . 0 , 

0 . 0 , 

0 . 0 , 

0 . 0 ); 


/* Special matrix of coefficients for time = 0 to 10 seconds */ 
/* added for start transients 6/26/91 by Mark O' Valentine: */ 


static float YNOM_START [22] [4] = ( 

0 . 0 , 

o.o, 

0 . 0 , 

0 . 0 , 

0 - 0 , 

0 . 0 , 

0 . 0 , 

0.0 

0.0 

0 . 0 , 

0 . 0 , 

0 . 0 , 

0 . 0 , 

0 . 0 , 

0 . 0 , 

0 . 0 , 

0 . 0 , 

0 . 0 , 

0 . 0 , 

0 . 0 , 

0 . 0 , 


0.0, 520.0420, 1777.250, -1176.211, 


520.0420, 

1777.250, 

-1176.211, 

-537.0478, 

2449.364, 

-1004.361 

-537.0478, 

2449.364 

, -1004.361 

0.0, 

0.0, 

o.o. 

o.o, 

o 

o 

o.o. 

o.o. 

o 

o 

o.o. 

o.o. 

o.o. 

o.o, 

93179.35, 

-124648.3, 

63325.11, 

93179.35, 

-124648.3, 

63325.11, 

o.o. 

0.0, 

0.0, 

o.o. 

0.0, 

0.0, 

4604.870, 

187.0750, 

554.5640, 

4604.870, 

187.0750, 

554.5640, 

4952.407, 

-1274.236, 

1205.108, 

4952.407, 

-1274.236, 

1205.108, 

3003.940, 

196.7420, 

o.o. 

3003.940, 

196.7420, 

o 

• 

o 

3006.000, 

0.0, 

© 

o 

3006.000, 

0.0, 

o.o. 

3006.000, 

0.0, 

o.o. 

3006.000, 

0.0, 

0.0>; 


/* Matrix of polynomial coefficients for model simulation; */ 
/* Mark's altered (improved) coefficients, circa 2/20/91: V 


static float COEF [22] [7] [4] = < 

6.75581E-02, 3.13551E-01, -6.33715E-01, 
-1.70377E-02, 1.42833E-02, -4.50418E-03, 
1.33900E-02, -3.66618E-03, -3.06278E-03, 
4.45270E-01 , -5.47878E-01, 3.84827E-01, 
7.81299E-02, -3.78262E-01 , 1.77979E-01, 


O.OOOOOE+OO, 

O.OOOOOE+OO, 

O.OOOOOE+OO, 

O.OOOOOE+OO, 

0.00000E+00, 
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2.48217E-03, 

-6.22494E-04, 

6.75581E-02, 

-1.70377E-02, 

1.33900E-02, 

4.45270E-01, 

7.81299E-02, 

2.48217E-03, 

-6.22494E-04, 

-1.79474E-01, 
-9.45672E-04, 
-3.15037E-02, 
3.09397E-02, 
-4.63672E-01 , 
-5.18391E-03, 
1.87290E-03, 

-1.79474E-01, 

-9.45672E-04, 

-3.15037E-02, 

3.09397E-02, 

-4.63672E-01, 

-5.18391E-03, 

1.87290E-03, 

-3.27388E+00, 

-2.72697E-02, 

4.14552E-02, 

1.65703E+00, 

-7.63844E-02, 

-1.04713E-02, 

-2.40354E-03, 

-3.27388E+00, 

-2.72697E-02, 

4.14552E-02, 

1.65703E+00, 

-7.63844E-02, 

-1.04713E-02, 

-2.40354E-03, 

1.33173E+00, 

6.04503E-04, 

1.30019E-02, 

2.20889E-02, 

8.40340E-01. 

-1.19188E-03, 

-9.63561E-04, 

1.33173E+00, 

6.04503E-04, 

1.30019E-02, 

2.20889E-02, 

8.40340E-01, 

-1.19188E-03, 

-9.63561E-04, 

-6.31739E-01 , 


-6.37037E-03, 

-3.36372E-04, 

3.13551E-01, 

1.42833E-02, 

-3.66618E-03, 

-5.47878E-01, 

-3.78262E-01, 

-6.37037E-03, 

-3.36372E-04, 

5.30960E+00, 
1.75486E-02, 
-4.35284E-02, 
-2.52552E-01 , 
1.84179E+00, 
1.99238E-02, 
3.37052E-03, 

5.30960E+00, 

1.75486E-02, 

-4.35284E-02, 

-2.52552E-01, 

1.84179E+00, 

1.99238E-02, 

3.37052E-03, 

8.50257E+00, 

6.61322E-02, 

-8.71833E-02, 

-4.27165E+00, 

2.17812E-01, 

2.78641E-02, 

6.10207E-03, 

8.50257E+00. 

6.61322E-02, 

-8.71833E-02, 

-4.27165E+00, 

2.17812E-01, 

2.78641E-02, 

6.10207E-03, 

-3.68307E+00, 

O.OOOOOE+OO, 

-2.74375E-02, 

-6.10283E-02, 

-2.36836E+00, 

2.29577E-03, 

2.13802E-03, 

-3.68307E+00, 

O.OOOOOE+OO, 

-2.74375E-02, 

-6.10283E-02, 

-2.36836E+00, 

2.29577E-03, 

2.13802E-03, 

1 . 12647E+00, 


2.89615E-03, 

3.61408E-04, 

-6.33715E-01, 
-4.50418E-03, 
-3.06278E-03, 
3.84827E-01 , 
1.77979E-01, 
2.89615E-03, 
3.61408E-04, 

-1. 78281 E+00, 
-1.05341E-02, 
4.25821E-02. 
O.OOOOOE+OO, 
-7.887&8E-01 , 
-8.76068E-03, 
-2.36554E-03, 

-1.78281 E+OO, 
-1.05341E-02, 
4.25821E-02, 
O.OOOOOE+OO, 
-7.88768E-01, 
-8.76068E-03, 
-2.36554E-03, 

-6.35848E+00, 

-4.57526E-02, 

5.34435E-02, 

2.88507E+00, 

-1.79975E-01, 

-2.02631E-02, 

-4.52982E-03, 

-6.35848E+00, 

-4.57526E-02, 

5.34435E-02, 

2.88507E+00, 

-1.79975E-01, 

-2.02631E-02, 

-4.52982E-03, 

3.04542E+00, 

O.OOOOOE+OO, 

O.OOOOOE+OO, 

O.OOOOOE+OO, 

1.88330E+00, 

O.OOOOOE+OO, 

O.OOOOOE+OO, 

3.04542E+00, 

O.OOOOOE+OO, 

O.OOOOOE+OO. 

O.OOOOOE+OO, 

1.88330E+00, 

O.OOOOOE+OO, 

O.OOOOOE+OO 

-8.22106E-01 


O.OOOOOE+OO, 

O.OOOOOE+OO, 

O.OOOOOE+OO, 

O.OOOOOE+OO, 

O.OOOOOE+OO, 

O.OOOOOE+OO, 

O.OOOOOE+OO, 

O.OOOOOE+OO, 

O.OOOOOE+OO, 

O.OOOOOE+OO, 

O.OOOOOE+OO, 

O.OOOOOE+OO, 

O.OOOOOE+OO, 

O.OOOOOE+OO, 

O.OOOOOE+OO, 

O.OOOOOE+OO, 

O.OOOOOE+OO, 

O.OOOOOE+OO, 

O.OOOOOE+OO, 

O.OOOOOE+OO, 

O.OOOOOE+OO, 

O.OOOOOE+OO, 

O.OOOOOE+OO, 

O.OOOOOE+OO, 

O.OOOOOE+OO, 

O.OOOOOE+OO, 

O.OOOOOE+OO, 

O.OOOOOE+OO, 

O.OOOOOE+OO, 

O.OOOOOE+OO, 

O.OOOOOE+OO, 

O.OOOOOE+OO, 

O.OOOOOE+OO, 

O.OOOOOE+OO, 

O.OOOOOE+OO, 

O.OOOOOE+OO, 

O.OOOOOE+OO, 

O.OOOOOE+OO, 

O.OOOOOE+OO, 

O.OOOOOE+OO, 

O.OOOOOE+OO, 

O.OOOOOE+OO, 

O.OOOOOE+OO, 

O.OOOOOE+OO, 

O.OOOOOE+OO, 

O.OOOOOE+OO, 

O.OOOOOE+OO, 

O.OOOOOE+OO, 

O.OOOOOE+OO, 

O.OOOOOE+OO, 

O.OOOOOE+OO, 

O.OOOOOE+OO. 
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-1 . 17272E-02 
1.19841E-03 
3.06835E-01 
-2.97835E-02 
* 1 .32908E-04 
-1.11418E-04 

-6.31739E-01 ( 
-1 . 17272E-02, 
1 . 19841E-03, 
3.06835E-01 , 
-2.97835E-02, 
-1.32908E-04, 
-1.11418E-04, 

-1.53861E-01, 
7. 03175E-05, 
-5.88139E-02, 
-4.07650E-03, 
1 . 13350E-01 ( 
-9.45481E-04, 
8.45892E-04, 

-1 .53861E-01 t 
7.03175E-05 ( 
-5.88139E-02, 
-4.07650E-03, 
1.13350E-01, 
-9.45481E-04, 
8.45892E-04, 

-1.06276E-02, 
-3.09662E-04, 
-1 .11283E-02, 
7.83933E-03, 
•2.95686E-02, 
•4.82056E-04, 
6.05011E-04, 

-1.06276E-02, 

-3.09662E-04, 

-1.11283E-02, 

7.83933E-03, 

-2.95686E-02, 

-4.82056E-04, 

6.05011E-04, 

-4.17196E-01, 
-3.58206E-03, 
1 .29300E-03, 
1.15425E-01, 
-2.13535E-02, 
-1 .19488E-04, 
-8.58883E-05, 

-4.17196E-01 , 
-3.58206E-03, 
1.29300E-03, 
1. 154256-01, 
-2.13535E-02, 


1.26742E-02 

-4.42708E-03 

-3.73406E-01 

8.02084E-02 

O.OOOOOE+OO 

3.62329E-Q4 

1 . 12647E+00 
1.26742E-02 
-4.42708E-03 
-3.73406E-01 
8. 02084 E- 02 
O.OOOOOE+OO, 
3.62329E-04, 

7.19883E-0V 
O.OOOOOE+OO, 
9.10559E-02, 
7.66015E-03, 
3.59439E-01 , 
4.14947E-03, 
O.OOOOOE+OO, 

7.19883E-01 , 
O.OOOOOE+OO, 
9.10559E-02, 
7.66015E-03, 
3.59439E-01 , 
4.14947E-03, 
O.OOOOOE+OO, 

4.20930E-02, 
-3.51086E-04, 
5.70Q60E-03, 
-2.56483E-04, 
2.10816E-01 4 
1.94540E-03, 
-1 .20095E-04, 

4.20930E-02, 
-3.51086E-04, 
5.70Q60E-03, 
-2.56483E-04, 
2. 10816E-01 ( 
1.94540E-03, 
-1.20095E-04, 

3.95345E-01 t 
2.19911E-03, 
-4.85527E-03 1 
-1.52808E-01, 
5.46972E-02, 
O.OOOOOE+OO, 
3.00620E-04, 

3.95345E-01, 
2.19911E-03, 
-4.85527E-03, 
•1 .52808E-01 g 
5.46972E-02, 


-5.26788E-03 

2.67544E-03 

2.76127E-01 

-4.00085E-02 

O.OOOOOE+OO 

-2.01927E-04 

-8.22106E-01 
-5.26788E-03 
2.67544E-03 
2.76127E-01 , 
-4.00085E-02, 
O.OOOOOE+OO, 
-2.01927E-04, 

-3.08034E-01 , 
O.OOOOOE+OO, 
-4.02095E-02, 
-5.73094E-03, 
-1 .50146E-01 , 
-1 .81101E-03, 
O.OOOOOE+OO, 

-3.08034E-01, 
O.OOOOOE+OO, 
-4. 020956-02, 
-5.73094E-03, 
-1.50146E-01, 
-1.81101E-03, 
O.OOOOOE+OO, 

O.OOOOOE+OO, 
2.48367E-04, 
O.OOOOOE+OO , 
1 .02111E-02, 
-7.71641E-02, 
-7.86290E-04, 
O.OOOOOE+OO , 

O.OOOOOE+OO , 
2.46367E-04, 
O.OOOOOE+OO , 
1.02111E-02g 
-7. 71641E-02, 
-7.86290E-04, 
O.OOOOOE+OO, 

-3.26333E-01 , 
-6.03413E-04, 
2.59230E-03, 
1 . 16439E-01 , 
*1 .49099E-02, 
O.OOOOOE+OO, 
-1.28705E-04, 


O.OOOOOE+iO, 

O.OOOOOE+OO, 

O.OOOOOE+OO, 

O.OOOOOE+OO, 

O.OOOOOE+OO, 

O.OOOOOE+OO, 

O.OOOOOE+OO, 

O.OOOOOE+OO, 

O.OOOOOE+OO, 

O.OOOOOE+OO, 

O.OOOOOE+OO, 

O.OOOOOE+OO, 

O.OOOOOE+OO, 

O.OOOOOE+OO, 

O.OOOOOE+OO, 

O.OOOOOE+OO, 

O.OOOOOE+OO, 

O.OOOOOE+OO, 

O.OOOOOE+OO, 

O.OOOOOE+OO, 

O.OOOOOE+OO, 

O.OOOOOE+OO, 

O.OOOOOE+OO, 

O.OOOOOE+OO, 

O.OOOOOE+OO, 

O.OOOOOE+OO, 

O.OOOOOE+OO, 

O.OOOOOE+OO, 

O.OOOOOE+OO, 

O.OOOOOE+OO, 

O.OOOOOE+OO, 

O.OOOOOE+OO, 

O.OOOOOE+OO, 

O.OOOOOE+OO, 

O.OOOOOE+OO, 

O.OOOOOE+OO, 

O.OOOOOE+OO, 

O.OOOOOE+OO, 

O.OOOOOE+OO, 

O.OOOOOE+OO, 

O.OOOOOE+OO, 

O.OOOOOE+OO, 

O.OOOOOE+OO, 

O.OOOOOE+OO, 

O.OOOOOE+OO, 

O.OOOOOE+OO, 

O.OOOOOE+OO, 

O.OOOOOE+OO, 


-3.26333E-01 

-6.03413E-04 

2.59230E-03 

1.16439E-01 

-1.49099E-02 


O.OOOOOE+OO, 
O.OOOOOE+OO, 
O.OOOOOE+OO, 
O.OOOOOE+OO, 
O.OOOOOE+OO, 
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0 . OOOOOE+OO f 0 . OOOOOE+OO , 
3.00620E-04, -1 .28705E-04, 


-1 .19488E-04, 
-8.58883E-05 , 

-2.66681E-02, 
-2.37199E-04, 
-3.22731E-04, 
1 .00202E-02 , 
1.17278E-03, 
2.65487E-06, 
2.85558E-05, 

-2.66681E-02, 
-2.37199E-04, 
-3.22731E-04, 
1.00202E-02, 
1.1 7278E - 03 , 
2.65487E-06, 
2.85558E-05, 


9.92401E-03, 

O.OOOOOE+OO, 

O.OOOOOE+OO, 

-1.39743E-02, 

-3.27962E-05, 

-1.72772E-04, 

-1.72770E-05, 

9.92401E-03, 
O.OOOOOE+OO, 
O.OOOOOE+OO, 
*1 .39743E-02, 
-3.27962E-05, 
-1.72772E-04, 
-1.72770E-05, 


-2.27310E-02, 

0.00000E+00, 

O.OOOOOE+OO, 

1.31286E-02, 

4.38514E-03, 

3.49690E-05, 

1.68002E-05, 

-2.27310E-02, 

O.OOOOOE+OO, 

O.OOOOOE+OO, 

1.31286E-02, 

4.38514E-03, 

3.49690E-05, 

1.68002E-05, 


O.OOOOOE+OO, 

O.OOOOOE+OO, 

O.OOOOOE+OO, 

O.OOOOOE+OO, 

O.OOOOOE+OO, 

O.OOOOOE+OO, 

O.OOOOOE+OO, 

O.OOOOOE+OO, 

O.OOOOOE+OO, 

O.OOOOOE+OO, 
O.OOOOOE+OO, 
O.OOOOOE+OO, 
O.OOOOOE+OO, 
O.OOOOOE+OO, 
O.OOOOOE+OO, 
O.OOOOOE+OO >; 


/*♦ 

/** 


MAIN BODY OF CALL_RSA() 


**/ 


/** (otherwise known as the control module for USUNO) **/ 



void call_RSA( input, input_tags, output, output_tags) 


struct inp 
struct inp_tags 
struct outp 
struct outp_tags 


♦input; 

*input_tags; 

♦output; 

* out put _ tags; 


< 


void do_nothing(); 

void perform_start_transient_model lingO; 

void clear_the_accufnulators_average_H_and_average_L(); 

void Calculate_the_weighted_si»n(); 

vo i d addas i ng l e_sample_to_f i rsti ni t_po i nt_average< ) ; 
void Calc_1st_avg_signature(); 

void Load_f i rst_signature_into__model_for_onej>oint_ini t(); 

void initialize_variables<); 
void Transfer_to_ADT(); 

indata = input; /* AWFUL KLUGE to get passed parameters to be GLOBAL. */ 

outdata = output; /* Sorry folks, but this program was inherited from a */ 
/* BASIC programmer. */ 


out_tags = output_tags; 

Trensfer_to_ADTO; /* Transfer data from parameter list to a local array */ 
/* for looping, indexing, processing (ADT) V 


/* Local variables are initialized separately in SAFO.C and/or in the V 
/* SAFO control module at Hunstville using extern void init(). V 

165 


/* This is done at the beginning of each new engine test. */ 

if (ADT £90] > 0.8) 
dojiothingO; 

if (very_f i rst_data _point == 1) 
i 

perf orm_s t ar t_t rans i entjnode U i ng( ) ; 
very_first_dat appoint = 0; 

> 

if (ADT [90] < time_that_we_begin_sampling) 

Ca l cul ate_the_we i ghted_sun( ) ; 

if ((ADTC90) >= time_that_we_begin_sampling) && (ADT [90] < end_of_sampl i ng_per i od ) ) 

< 

add_a_s i ng l esampl e_to_f i rst_ini t_poi nt_average( ) ; 

$ample_count = samplecount ♦ 1; 

Ca l cu l a te_the_we i ghted_sun( ) ; 

> 

if ( (ADT (90] >= endofsampl ing_period) U ( donees amp l ing « 0)) 

< 

Ca l c_t st_avg_s i gnature( ) ; 

Load_f i rs t_s i gna t ure_ i n t ojnode l Jt or_one _po i nt_ i n i t ( ) ; 
done_sampl ing = 1; 

> 

if (done_sampl ing 1) 

Ca l cu l a t e_ t he_we i gh t ed_s un( ) ; 

Pc_old « ADT £34] ; /* Preserve history data for later calcs. */ 

> 


/** **/ 

/** RSA SUBROUTINES **/ 

/** **/ 

/** A subroutine called by one of these belongs to the "E M calling level. **/ 
/** **/ 


void do_nothing() 

C 

/* used for debugging under the PC environment Turbo C */ 

> 


int init() 

< 

int i. a. b, c; 


/* Clear out the weighed sum calculation variables: */ 
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for (i = 0; i < 11; i++) 
npdli] = 0.0; 
for (i = 0; i < 22; i 

estli] = 0.0; 
ynomt [il = 0.0; 

> 


/* Clear out the accumulators where the samples are stored: */ 

for (i=0; i<max_adt; i++) 
average_H[i]=0; 
for (i=0; i<max_adt; i++) 
average_l [i] =0; 

wfax_tot = 0.0; /* clear out the backpressure suns */ 

woax_tot = 0.0; 

for < i =0; l <22; i++) 

deltay L[i] = 0.0; /* clear out the 5-second delta-y sun */ 

/* Initialize the filtering array (created by Ed Nemeth 9-3-92) V 

for (a = 0; a < 22; a++) 
for (b = 0; b < 10; b++) 
f i lt_nd[al [bl = 0.0; 


/* initialize working variables: */ 

end_of_sampl ing_per i od = time_to_wai t+ tiniest ha t_we_begin_samp ling; 
sample_count = 0; 
very_f i rst_data_point = 1; 
done_sarcpling = 0; 

return<1); 

> 

void Transfer_to_ADT( ) 

ADT [0] = indata->adt0; 

ADT 11] = indata->adt1; 

ADT [2] = indata->adt2; 

ADT [3] = indata->adt3; 

ADT [4] = indata->adt4; 

ADT [5] = indata->adt5; 

ADT [6] * indata->adt6; 

ADT [7] = indata->adt7; 

ADT [81 = indata->adt8; 

ADT [9] = indata->adt9; 

ADT [101 = indata->adt10; 

ADT [11] = indata->adt11; 

ADT [121 = indata->adt12; 

ADT [13] = indata->adt13; 

ADT [14] = indata->adt14; 

ADT [15] = indata->adt15; 

ADT [16] = indata->adt16; 

ADT [17] = indata->adt17; 

ADT [18] = indata->adt18; 

ADT [19] = indata->adt19; 

ADT [20] = indata->adt20; 

ADT [21] = indata->adt21; 
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ADT [22] = indata->adt22 
ADT [34] = indata- >adt34 
ADT [38] = indata- >adt38 
ADT [39] = indata->adt39 
ADT [40] = indata- >adt40 
ADT [41] = indata->adt41 
ADT [42] = indata- >adt42 
ADT [43] = indata- >adt43 
ADT [47] = indata- >adt47 
ADT [48] = indata- >adt48 
ADT [49] = indata- >adt49 
ADT [50] = indata->adt50 
ADT [51] = indata->adt51 
ADT [52] = indata->adt52 
ADT [53] = indata->adt53 
ADT [54] = indata->adt54 
ADT [55] = indata->adt55 
ADT [56] = indata->adt56 
ADT [57] = indata->adt57 
ADT [90] = indata- >adt90 


void perform_start_transientjnodel l ing< ) 

/* Models startup portion of test from 0.2 to 10.1 seconds by using */ 

/* modified YN0M array with special values in it derived from actual */ 

/* data values. See variable YN0H_START in include file NH0D3.H. V 

/* Mark O'Valentine - formulas; Robert Van den Meuvel - coding 7/10/91 V 

int i; 
float PL; 

PL = ((ADT [19]+ADT [20] )/2)/3006; 
for ( i=0; i <22; ?♦+> 

< 

YNOH_START [i] [0] = ADT [i+1] - 

(YNOM_START[i] [1]*PL) - 
(YN0M_START [i] [2]*(PL*PL)) - 
(YNOM__START[i] [3]*(Pl*PL*PL)); 

v[i] = 1.0; /* The v-array tells whether a measurement is valid. V 

) 

/* be sure not to set HPFT speeds to the uncleared, bogus values */ 

/* (about 40,000 rpm) from the pre-test sensor checkout */ 

/* (correct speeds should not be > 4000 rpm at 0.2 seconds V 

if (ADT [9] > 4000) 

C 

YM0M_START [8] [0] = 456.0 - 
(YN0H_START[8J C1]*PL) - 
(YN0H_START[8] C2]*(PL*PL)) - 
(YNGH_START [8] [3]*(PL*PL*PL)); 

> 

if (ADT [10] > 4000) 
i 

YMOM_START [9] [0] = 456.0 - 
(YN0M_$TART [9] [1]*PL) - 
(YM0M_START[9] [2]*(PL*PL)> - 


168 



(YNOM_START [9] £3]*(PL*PL*PL)); 

> 

> 


vo i d add_a_s i ng l e_samp l e_t o_f i r s t_i ni t_po i nt_a verage( ) 

/* Takes one raw data point and adds it to a naming average which will V 
/* be used later to produce an engine signature at the current pwr level. */ 

i 

int i; 

float wfaxO, woax(); 

for (i =0; i < max_adt; i++) 

average_H£i] = average_H [i]+ ADT M3; 

wfax_tot = wfaxtot + wfaxO; /* these are for same purpose V 

woax tot = woax_tot ♦ woax(); 


void Calc_1st_avg_signature() 

/* M sampte_count" counts the number of data points sampled so far. */ 
/* Ue average the last 5 seconds of raw data to get an average for V 
/* each parameter ( average Jifx3 ). Each average is a "signature” for V 
/* the engine at the current power level. V 

( 

int i; 

int item[7] = (0,38,47,41,50,0,0}; 
float PL, PL2, PL3; 

for <i=0; i<max_adt; i«“0 

average_H[i] = average_HCi] / sample_count; 

PL * average_H £343/3006; 

PL2 = PL * PL; 

PL3 = PL * PL2; 

/* set baseline indep. params to their 1st init. point average values */ 


XNOME03E0] = 0.0; 
for (i=1; i<7; i++) 

( 

ED_N0M£i3 = XNOMEHE0] ♦ 

XNOM [ i ] [1] * PL + 

XNOM [i] [2] * PL2 + 

XN0M£i3 £3] * PL3; 

> 

for (i=1; i<5; i++) 

( 

XNOM £ i ] £03 = (XNOM Ei] £03 * average_H£i tem£i33 )/ED_N0M£i3 ; 
XNOM [i] £13 = (XN0M£i3E13 * average_H£i tem£i)3 )/ED_N0M£i3 ; 
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> 

i 


1 


XNOHCi] [21 = (XNOHCi] t<'3 * averagejl [itemCi]] )/ED_NOHCi] ; 
XNOHCi] C3] * (XNOHCi] [3] * average_H (i temCi] ] )/ED_NOM (i ) ; 


XNOHCi] CO] 
XNOHCi] Cl] 
XNOHCi] C2J 
XNOHCi] C3] 

XNOHCi] CO] 
XNOHCi] Cl] 
XNOHCi] C2] 
XNOHCi] C3] 


(XNOHCi] CO] * (wfax_tot/sample_count) )/ED_NOHCi] ; 
(XNOHCi] Cl] * (wfax_tot/sample_count))/ED_NOHCi]; 
(XNOHCi] [2] * (wfax_tot/safflple_count))/ED_NOH(i] ; 
(XNOHCi] [3] * (wfax_tot/sanple_count))/ED_NOH(i] ; 

(XNOHCi] CO] * (woax_tot/sample_count))/ED_NOHCi]; 
(XNOHCi] Cl] * (woax_tot/sample_count))/ED_N(JHCiJ; 
(XNOHCi] [2] * (woax_tot/sample_count ))/ED_NOH [i] ; 
(XNOHCi] [3] * (woax_tot/sample_count))/ED_NOHCi J ; 


/* debug 7-24-92 for Ed Nemeth: Failure investigation. Rewrite of XNOH 
calculation to custom specs for debugging. RVDH 

for (i=1; i<5; i++) 

XNOHCi) CO] = average_H C i tern C i ] ] - XNOHCi] C1]*PL - 
XNOHCi] C2] * PL2 - XNOHCi) C3] * PL3; 

XNOH [5] CO] = wfax_tot / sample count • XNOH C5] CD * PL - 
XNOHCS] C2] * PL2 - XNOHCS] C3] * PL3; 

XNOH [6] [0] = woax.tot / sample_count - XN0H[6]CD * PL - 
XN0HC6] [2] * PL2 - XN0HC6] C3] * PL3; 

*/ 


void Load_f i rst_s i gnature_i nto_mode l_f or_one_point_i ni t ( ) 

C 

float PL1, PL2, PL3; 
int i; 

PL1 = average_HC34]/3006; /* define the current power level t powers of V 

PL2 = PL1 * PL1; 

PL3 = PL1 * PL1 * PL1; 

/* LORD THE 104X SIGNATURE INTO THE YNOH HATRIX: */ 
for (i=0; i < 22; i++) 

C 

/• YNOH [i] [0] IS THE 1ST COLUNN OF A 22x4 HATRIX. THE 1ST COLUHN */ 

/* DEFINES THE CONSTANT TERNS IN A SET OF 22 3RD DEGREE POLYNOH1ALS. */ 

/* VH = "V-High Average, •• a variable which represents the selected, */ 

/* UNADJUSTED V-High‘s for 1-point equation calculations */ 


VHCi] = average_HCi+1] ; 

YNOHCil CO] = VHCi] - PL1*YN0HCi) CD - PL2*YNOHCi] C2] - PL3*YNOHCi] C3] ; 


> 

) 
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void Calculate_the_weighted_sum( ) 
i 

float USUMO; 

float w_sum = -100; 

w sum * WSUNO; /* Redlining sum (weighted) is calculated here 

f* TRANSFER ALL PARAMETERS OF INTEREST TO THE OUTPUT ARRAYS: V 

outdata->adtO - ADT [0] ; 
outdata- >adt1 = ADT [13; 
outdata->adt2 = ADT [23? 
outdata->adt3 = ADT [33 ; 
outdata->adt4 = ADT [43 ; 
outdata- >adt5 = ADT [53; 
outdata->adt6 = ADT [6]; 
outdata->adt7 = ADT [73; 
outdata- >adt8 = ADT [83; 
outdata- >adt9 = ADT [93; 
outdata- >adt 10 = ADT [103; 
outdata->adt11 = ADT [113; 
outdata->adt12 = ADT [123; 
outdata- >adt 13 = ADT [133 ; 
outdata->adt14 = ADT [143; 
outdata->adt15 = ADT [153; 
outdata->adt16 = ADT [163; 
outdata- >adt 17 = ADT [173; 
outdata->adt18 = ADT [183; 
outdata->adt19 = ADT [193; 
outdata ->adt20 = ADT [203; 
outdata- >adt21 = ADT [21 3 ; 
outdata->adt22 = ADT [223; 
outdata- >adt34 = ADT [343; 
outdata- >adt38 = ADT [383; 
outdata- >adt39 = ADT [393; 
outdata- >adt40 = ADT [403; 
outdata->adt41 = ADT [413; 
outdata->adt42 = ADT [423; 
outdata- >adt43 = ADT [433; 
outdata- >adt47 = ADT [473; 
outdata- >adt48 = ADT [483; 
outdata- >adt49 = ADT [493; 
outdata- >adt50 = ADT [503; 
outdata- >adt51 * ADT [513; 
outdata->adt52 = ADT [523; 
outdata- >adt53 = ADT [533; 
outdata ->adt54 = ADT [543; 
outdata- >adt55 = ADT [553; 
outdata- >adt56 = ADT [563; 
outdata->adt57 = ADT [573; 
outdata- >adt 90 = ADT [903; 
outdata- »est0 * est [03 ; 
outdata- >est1 = est [13; 
outdata- >est2 = est [23; 
outdata->est3 = est [33; 
outdata ->est 4 = est [43; 



outdata->est5 = est [5] ; 
outdata->est6 * est [6]; 
outdata->est7 = cst (7]; 
outdata->est6 = estC8}; 
outdata->est9 = est[9]; 
outdata->est10 * est [10]; 
outdata->est11 = est [11]; 
outdata->est12 = est [12]; 
outdata->est13 = est [13]; 
outdata->est14 = est [14]; 
outdata->est15 = est[15]; 
outdata->est16 = est [16]; 
outdata->est17 = est [17]; 
outdata->est18 = est [18]; 
outdata->est19 = est [19]; 
outdata->est20 = est [20] ; 
outdata->est21 = est [21] ; 
outdata->wsun = w_sun; 
outdata->shutdown = 0; 

f* SET FLAGS TO HEX ZERO SINCE WE DON'T USE THIS FEATURE: •/ 

out_tags->adtO * OxOOOOOOOOL; 
out_tags->adt1 = OxOOOOOOOOL; 
out_tags->adt2 = OxOOOOOOOOL; 
out_tags->adt3 = OxOOOOOOOOL; 
out_tags->adt4 = OxOOOOOOOOL; 
out_tags->adt5 = OxOOOOOOOOL; 
out_tags->adt6 = OxOOOOOOOOL; 
out_tags->adt7 = OxOOOOOOOOL; 
out_tags->adt8 = OxOOOOOOOOL ; 
out_tags->adt9 = OxOOOOOOOOL; 
out_tags->adt10 = OxOOOOOOOOL 
out_tags->adt11 = OxOOOOOOOOL 
out_tags->adt12 = OxOOOOOOOOL 
out_tags->adt13 = OxOOOOOOOOL 
out_tags->adt14 * OxOOOOOOOOL 
out_tags->adt15 = OxOOOOOOOOL 
out_tags->adt16 = OxOOOOOOOOL 
out_tags->adt17 = OxOOOOOOOOL 
out_tags->adt18 = OxOOOOOOOOL 
out_tags->adt19 = OxOOOOOOOOL 
out_tags->adt20 * OxOOOOOOOOL 
out_tags->adt21 = OxOOOOOOOOL 
out_tags->adt22 = OxOOOOOOOOL 
out_tags->adt34 = OxOOOOOOOOL 
out_tags->adt38 = OxOOOOOOOOL 
out_tags->adt39 = OxOOOOOOOOL 
out_tags->adt40 = OxOOOOOOOOL 
out_tags->adt41 = OxOOOOOOOOL 
out_tags->adt42 = OxOOOOOOOOL 
out_tags->adt43 = OxOOOOOOOOL 
out_tags->adt47 = OxOOOOOOOOL 
out_tags->adt48 = OxOOOOOOOOL 
out_tags->adt49 * OxOOOOOOOOL 
out_tags->adt50 = OxOOOOOOOOL 
out_tags->adt51 = OxOOOOOOOOL 
out_tags->adt52 = OxOOOOOOOOL 
outtags- >adt53 = OxOOOOOOOOL 
out_tags->adt54 = OxOOOOOOOOL 
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out_tags->adt55 = 0 x 000000 00 L; 
out_tags - >adt 56 = OxOOOOOOOOL; 
out_tags->adt57 = OxOOOOOOOOL; 
out_tags*>adt90 = OxOOOOOOOOL; 
out_tags->estO = OxOOOOOOOOL; 
out_tags->est1 = OxOOOOOOOOL; 
out_tags->est2 = OxOOOOOOOOL; 
out tags->est3 = OxOOOOOOOOL; 
out_tags->est4 = OxOOOOOOOOL; 
out_tags->est5 = OxOOOOOOOOL; 
out_tags->est6 = OxOOOOOOOOL; 
out_tags->est7 = OxOOOOOOOOL; 
out tags->est8 = OxOOOOOOOOL; 
out_tags->est9 = OxOOOOOOOOL; 
out_tags->est10 = OxOOOOOOOOL; 
out_tags->est11 = OxOOOOOOOOL; 
out tags->est12 = OxOOOOOOOOL; 
out_tags->est13 = OxOOOOOOOOL; 
out tags->est14 - OxOOOOOOOOL; 
out_tags->est15 = OxOOOOOOOOL; 
out_tags->est16 = OxOOOOOOOOL; 
out — tags->cst17 = OxOOOOOOOOL; 
out_tags->est18 = OxOOOOOOOOL; 
out_tags->est19 = OxOOOOOOOOL; 
out tags->est20 = OxOOOOOOOOL; 
out_tags->est21 = OxOOOOOOOOL; 
out_tags->wsun = OxOOOOOOOOL; 
out tags->shutdown = OxOOOOOOOOL; 


> 


'**/ 

*+/ 


/+* WEIGHTED SUH SUBROUTINE *V 

/** 

/** (calculates redlining sum of deviations between models and reality) *V 

**/ 

/** ... * v 
/** SUMMARY: This module is essentially one big list of equations which **/ 

/** get executed when the value of the redline sun is desired. **/ 

/** The module is cownented into modules (A, B, C, etc.) which **/ 

/*+ each represent some significant sub- process. **/ 



float WSUMO 
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/* VARIABLES USED FOR GENERAL WSUH PROCESSING */ 


< 


int alstef fensindex, alkey; 

float PL; 

float sun; 

int c,m,n,e,b,i, j; 

float 81, B2 f rl; 

int al, a2; 

/* Ed Nemeth: 9-3-92: variables added for filtering of nd calculation V 

float filt_sum = 0; 

int nsamp * 0; 

int nuraber_of_samples = 0; 

/* continuing with old WSUH definitions: V 

float AMIN = 0.2; /* not defined yet V 

float dmax = 0; 

float dl s 0, d2 = 0; 

float d [223 ; 

float vmaxMlj; 

int jq, kq, indx[3 ] ; f* added for Al Stef fan's modifications of MODULE C V 

float adiff[3], asum[3], aminimun = 3e+38; 

float nonredundant_pi delimit *11; /* Added 2/26/92 */ 


/* VARIABLES FOR WSUH WEIGHTING ( M W CALCULATION BELOW) added 2/92 */ 

float sig[22]; 
float sig_counter =0; 
float k(22] = < 

1, 1, 1.04, 1.13, 1.27, 1.49, 1.82, 2.31, 3.05, 4.20, 6.05, 8.91, 
13.91, 22.40, 37.76, 65.54, 120.15, 226.57, 440.56, 881.11, 1982.50, 
3965.00 >7 


float wfaxO; /* USED TO CALC (FUEL) REPRESSURIZATION FLOWS */ 
float woaxO; /* DITTO (LOX) V 


/* NIRD CODE CHANGES (01/16/92) 

we need to wait longer after power level transients are over 
before we begin computing PL based on commanded Pc again. 

if time < 5.9 sec, then compute PL using measured Pc: */ 


if (ADT [90] < time_main_stage_i screeched) 
i 

PL = ( (ADT [19] +ADT [20] )/2)/3006; 
transientcount = 0; 
transientf lag =1; 


/* else if time >= 5.9 sec, then check whether ADT [34] = Pc_old */ 
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else 

< 

if ( ADT [34] == Pc_old > 

< 


/+ if transient flag = 1 here, then it's the first sample after a 
transient, this code calculates PL using commanded Pc during 
steady-state engine operation. PL is calculated using measured 
Pc during transients, this code waits 4 samples after a transient 
is over before PL is computed using commanded Pc again. */ 

if ( transient_f lag == 1 ) 

< 

transient_count = 0; 

PL = ((ADT [19] +ADT 120] )/2)/3006; 

> 

/* else if transient^ lag = 0, then determine whether to compute */ 
/* PL using commanded Pc or measured Pc V 

else 

transient_count = t rans i ent_count ♦ 1; 

/* if transient_count < 4 then compute PL using measured Pc V 
if ( transient_count < 4 > 

PL = ( (ADT (193 +ADT [20] )/2)/3006; 

> 

/* if transient_count >= 4 then compute PL using commanded Pc V 
/* (this method is used during the majority of most tests) V 

else 

i 

PL = ADT [34] / 3006; 

> 

> 

/* since ADT [34) = Pc_old here, then re-set transient^ lag to 0 */ 
transient^ lag = 0; 

> 


/* else if ADT [34] <> Pcold then compute PL using measured Pc V 

else 

< 

PL = ((ADT (19) +ADT 120] )/2)/3006; 
t rans ient_f lag - 1; 
transient_count = 0; 

> 


/* M VAl M EQUATION CALCULATION: NODULE C V 



/* Validate six independent parameter measurements by comparison of V 

/* three redundant values V 
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i* Rewritten for speed and brevity by Al Steffens, Jr. 9-4-91 up to */ 
/* the statement, M val[n] = r; M */ 

for (n * 0; n < 6; n++) 

for (j<p0;jq<3; jq++) t 

/*** Get the ADT indices -- jq = 1st index, kq * off- index ***/ 
kq = (jq+1) X 3; /* permute kq */ 

indxtjq] * crossref [nit jq] ; /* ADT indices */ 

indxCkq] = cross ref [n] Ckq] ; /* ADT off-indices */ 

/*** Get ADT differences (and sums) ***/ 
adiffCjq] = fabs(ADT [indx[jq]] - ADT [indx[kq]] ); 
asumCjq] = (ADT [indxljq] ] ♦ ADT [ i ndx [kq] ] )/2; 

/*** Find smallest difference ***/ 
if (adiff[jq] < ami ni nun) { 
am ini nun = adiff[jq]; 
valCn] = asumljq]; 

> 

> /* for jq */ 

/*** Test if difference is within limit ***/ 
if (aminimm > PCHK[n]) 
valtn] = 0; 

> /* for n */ 

if (ADT [90] < timejnain_stage_is_reached) 

< 

val [6] = 0; 

> 

else 

val [6] = wfaxO; /* Call fuel repressurization flow calculation */ 

if (ADT [90] < time_main_stage_is_reached) 
i 

val [73 * 0; 

> 

else 

val [7] = woaxO; /* Call LOX/HEX repressurization flow calculation */ 

/* "INDVAL" EQUATION CALCULATION: */ 


/* */ 

/* Calculates mixture ratio (not yet implemented) and shifts indices/ values*/ 
/* to a new array (indval[73). The values in this new array are used later*/ 
/* in the “DELTA-X" calculation */ 

for (n * 0; n < 7; n++) 
i 

if (n == 0) 

indvaUn] = 0.0; /* mixture ratio */ 

> 

else 

if <(n == 5) || (n *= 6)) 

< 


/* repress */ 
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/* all others */ 


indval [n] = valtrHl]; 

> 

else 

indval [n] = val Cn] ; 

> 


/* ■•DEPNOH" EQUATION CALCULATION: HOOULE D */ 

/* v 

/* Calculate the nominal dependent values using PL in a polynomial V 
/* USE SPECIAL MATRIX YNOM_START (stands for YNOM- star tup- transient) */ 

/* THE TIME IS LESS THAN 5 SECONDS INTO THE TEST. MOOELS STARTUP */ 

/* SEPARATELY TO PROVIDE SMOOTHER WEIGHTED SUM. 6/27/91 RVdH */ 


if (ADT [90] < 10.0) 

L 

for (m = 0; m < 22; m*+) 

C 

depnomtmj = YN0M_$TART (m3 CO] + (YNOM_START Cm] [1]*PL) ♦ (YNOM_START [m3 [23 
* (PL * PL)) + (YNOM_START[m3 [3] * (PL * PL * PL)); 

> 

> 

else 

/* THE REGULAR MATRIX IS USED FOR THE REST OF THE TEST: V 

L 

for (m = 0; m < 22; nw-O 
( 

depnomlnO = YNOMIm][0] ♦ (YNOM M C1]*PL) + (YN0M[m][2] * (PL * PL)) ♦ (YNOM [m3 C3) * (PL * PL * PL)) 

> 

> 

for (m = 0; n < 22; m*+) 

< 

dep_star[m] = ADT [m+13/(1+deltay[m3 ); 

> 

/* "INDNOM" EQUATION CALCULATION: V 


/* Calculate the mominal independent values using PL in a polynomial V 

for (n =0; n < 7; n++) 

L 

indnomtn] = XNOMtnHO] + (XNOM[n3 tl) * PL) ♦ (XNOM[n][2l * (PL * PL)) 

♦ (XNOM [n] [3] * (PL * PL * PL)); 

> 

/* M DELTA"X M EQUATION CALCULATION: V 

/* v 

/* Calculate the X deltas between nominal and measured independent */ 

/* variable values */ 

for (n = 0; n < 7; n++) 

C 

B1 = indval [rO; /* indval array calculated in Module C V 

B2 = indnom[n3; 

if ((B1 == 0.0) || (B2 == 0.0)) 
rl = 0.0; 
else 
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rl = <B1 - B2)/B2; 

> 

del tax Cn] = rl; 

> 

/* •■INFLU" EQUATION CALCULATION: */ 

/* V 

/* Uses 30 influence coefficient array to derive estimates of dependent */ 
/* variables */ 

/* SET THE WHOLE COEFF ARRAY TO ZERO IF V 

/* THE TIME IS LESS THAN 5 SECONDS INTO THE TEST. HOOELS STARTUP */ 

/* SEPARATELY TO PROVIDE SMOOTHER WEIGHTED SUM. 6/27/91 RVdH */ 

if (ADT [90] < 10.0) 

/* IN THE BEGINNING, ZERO OUT THE COEFF ARRAY FOR STARTUP MODELLING: V 

< 

for (m = 0; m < 22; m++) 

for (n * 0; n < 7; m-O 

< 

influCmHn] = 0; 

> 

> 


> 

else /* FOR THE REST OF THE TEST, USE THE ARRAY AS NORMAL: */ 

i 

for (m = 0; m < 22; m++) 

< 

for (n = 0; n < 7; n++) 

< 

influCmHn] = COEFCmHrOCO] ♦ (COEFCmHnHI] * PL) ♦ 

(COEFCmJ Cn] [2] * (PL * PL)) + (COEFCm] Cn] C33 * (PL * PL * PL)); 

> 

> 

> 


/* M DELTA-Y M EQUATION CALCULATION: */ 

/* V 

/* Calculate the total influence of independent parameters on each */ 
/* dependent parameter V 


for (m = 0; m <22; m++) 
i 

sum = 0; 

for (n = 0; n < 7; n+*) 

sun = sum ♦ (influCmHn] * deltaxln]); 
de It ay Cm] = sum; 

> 


/* M EST" EQUATION CALCULATION: */ 

/* V 

/* Calculate the final estimate of the nominal dependent values V 


for (m = 0; m < 22; 

estCm] = depnomCm] * (1 + de It ay Cm] ); 


/* MOOULE E 

/* "D" EQUATION CALCULATION: Defined later, defaults used for now. 


*/ 

V 


/* No code here at this time (D'Val) 


10/30/91 RVDH V 


/* "V" EQUATION CALCULATION: HOOULE F */ 

*/ 


/* Initialize the valid "flag" array V 

/* if via] = 1 the sensor is valid and used in WSUM V 

/* if y[a] = 0.5 the sensor is temporarily invalid & used later V 


/* if v(a) = 0 the sensor is permanently invalid & never used again */ 

if (ADT [90) < 10.0) C 
for (a = 0; a < 22; a++) 

C 

d[a] = d_st la] ; 

if (a < 11) vmaxta) = vmax_st [a] ; 

if (vta) > 0) 

< 

vta] = 1.0; 

if (dial < 0) via) = 0.5; 
if (ADTta+1) == 0.0) vta) = 0.0; 

> 

else vta) = 0.0; 

> 

/* don't use HPFT speeds if their values are actually the uncleared V 
/* bogus values (-40,000 rpm) from the pre-test sensor checkout */ 

/* (correct speeds should not be > 4,000 rpm at 0.7 seconds */ 

if ((ADT [90] < 0.7) && (ADT [9] > 4000)) vt8) = 0.5; 
if ((ADT [90] < 0.7) U (ADT (101 > 4000)) v(9) = 0.5; 

/* don't use H.G. lnj Pressure (PID 24) if it's less than MCC Pc V 
/* after 1.5 seconds 

if ( (ADT C90] > 1.3) && (ADT [17] < PL*3006)) v(16) = 0.0; 

> 

else 

< 

for (a = 0; a < 22; a++) 

< 

d[a] = d ins [a] ; 

if (a < 11) vmax la) = vmaxjns [a] ; 


if (v(a) > 0) 

< 

vta) = 1.0; 

if (dta) < 0) vta) = 0.5; 
if (ADTta+1) == 0.0) vta) = 0.0; 

> 

else vta) = 0.0; 

> 

> 


/* Mark D'Val changes to ND calculation as of 9/10/91: 

*************************************************** 


V 

V 


/★ 41 MD ** EQUATION CALCULATION: 
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/* V 

/* Calculate normalized measurement deviations */ 

for (a = 0; a < 22; a++> 

if (via] < 1) 
edndla] = 0.0; 
else 

ednd la] = (ADTla+1] - est [a] )/d[a] ; 

> 

/* Added changes 9-3-92 by Ed Nemeth to smooth noise (filtering) */ 
number_of_samples = 5; 

n_samp = ntrnber_of _sampl es - 1; /* for C numbering system; starts at 0 V 

/* update filter array with latest data. Number_of_samples represents the 
number of buffer cells, each one of which represents a single sample 
in time, progressively back in time each data point observed. V 

for (a = 0; a < 22; a++) 

for (b = 0; b < nsamp; b*+) 

f i l t_nd [a] [b] = f i lt_r>d[a] [b+1] ; 

> 

for (a = 0; a < 22; a++) 

C 

f i l t_nd [a] Insamp] = ednd [a] ; 

> 

/* This double loop is to calc the average of Cnunber_of_samples] 
buffer cells as mentioned above. This represents the average of 
the last (number_of_samples] number of data points observed. */ 

for (a * 0; a < 22; e++> 

C 

f i lt_sum = 0; 

for (b * 0; b < nuraber_of _samp l es ; bM-) 

C 

filt_sum = filt_sum + filt_ndta] lb); 

> 

nd(a] = f i lt_sun/nunber_of_samples; /* nd is assigned the filtered value */ 


/* M VND" EQUATION CALCULATION: V 

/* V 


/* Calculate validated, normalized deviations of each measurement (0 - 21) */ 


for (a = 0; a < 11; a++) 
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al = 2.0 * a; 
a2 = al ♦ 1; 


/* accounts for dual sensors of some measurements */ 


/* cknax is the maximum standard deviation allowed between “good 11 V 
/* measured and estimated values V 

cknax = 3.0; /* set cknax to 3 until a better value is determined */ 

/* if both of the redundant parameters are available then check V 
/* whether one or the other should be disabled (never disable both V 
/* (if both redundant params indicate a problem use both in USUM */ 

if C(v[a1] == 1) M (v[a2) == D) 

< 

/+ if pi and P2 indicate similar trends (i.e., their measured values V 
/* are similar direction and magnitude from their estimated values), V 
/* then use both of them in the USUM calculation: 

if <fabs(nd[a1) * ndta23) < vmax[a]*dmax) 

( 

vial] = 1; 
v[a23 = 1; 

> 

else 

/* if Pi checks out but P2 does not, then disable P2 */ 

if ((fabs(nd[a1] ) < cknax) && ( fabs(nd[a2] ) > cknax)) 
v[a2) = 0; 
else 

/* if P2 checks out but PI does not, then disable PI V 
if ((fabs(nd[a13) > cknax) && (fabs(nd£a2) ) < cknax)) 
v[a1) = 0; 

> 

vndtal] = ndCal] * v[al3; 
vnd[a2) = nd[a23 * vta23; 


> 


/* validateO; */ 

/ft n HPD ii EQUATION CALCULATION: M00ULE G V 


/* calculate normalized parameter deviations for final summation of wsum V 

for (a=0; a < 11; a**) 

L 

al = 2.0 * a; 
a2 = al + 1; 

/* divide by the number of valid parameters, dl = 1 or 2 */ 

/* (set dl = 1 if vial] ♦ v[a2] <2 to avoid a divide by 0 or 1.5 */ 

dl = v[a13 ♦ vta23; 
if (dl < 2) dl = 1.0; 

npdta] = (fabs(vnd[a13) ♦ fabs(vnd[a23 ))/d1; 

> 



/* "W" EQUATION CALCULATION: */ 

/* v 

/* Define a weight for each normalized parameter deviation +/ 


sum = 0; /* set wsun to zero initially before calculating it */ 

/+ INITIALIZE THE SIG(] TABLE FOR WSUH WEIGHTING: (Added 2/92 RVdH ) : */ 

sigtO] = 0.0; 
sigCI] = 0.1; 
sig_counter =0.1; 
for ( i=2; i<22; i++) 

sig_counter = sig_counter + 0.2; 
sigti] - sig_counter; 

> 


/* CALCULATE THE WEIGHTS FOR WSUH BASED UPON A LINEAR FIT AIMED AT 
INCORPORATING NORMAL DISTRIBUTIONS RVdH/D' Valentine 2/92. */ 


/* 

for ( j=0; j<11; j++) 

< 

w[j] = 4000.00; 
for (i=0; i<21; i++) 

if ((npdtj) >= sigti)) && (npdtj) < sigti+1))) 

< 

wtj] = kti) + ( ((kti+1) - kti)>/0.2> * (npdtj) - sig(ij) ); 

npdtj) = 0; 

> 

> 

*/ 

/* CALCULATE THE WEIGHTED SUM OF DEVIATIONS FROM NORMAL MODEL VALUES: */ 

for (c=0; c<11; C++) 

< 

al = 2*c; 

•2 * «1 ♦ 1; 

wtc] = 1 + (npdtc]*npdtc]*npdlc)); 

if (((vtal) < 1) || (vta2) < 1)) U (wtc) > nonredundant_pid_ limit)) 
wtc] = nonredundant_pi delimit; 

sun = sun + fabs(wtc)); 

> 

return(sum); 

> /* exit WSUMO */ 


r**************** j 


LEVEL "F" SUBROUTINES 


/** 

/** 

/** 

/** 

/** 

/**' 


These are the bottom-most routines in the calling tree. 


**/ 

**/ 

**/ 

**/ 

**/ 

'**! 


float woax() 

( 

/* Calculate oxidizer backpressure. Mechanics unknown and assumed to work. V 
/* Robert Van den Heuvel (RSA programmer) 11/26/91 */ 

float hxpiv, hxtlv, hxdpv, hxc, betav, hxgamv; 
float gfl, gf2, gf3, gf4, hxy; 

float hxdiv, hxdtv, hxprat, hxb4, hxkc, hxtxf, hxrov; 
float pi,r, rl, r3; 
float sun; 

pi = 3.14159; 

hxpiv = ADT[56] + 14.7; 
hxtlv = ADTC57] ; 
hxdpv = ADTT58); 

/* Set Uoax parameters to reasonable values if the measured parameters */ 

/* are missing from the incoming data stream: (12/5/91 change -- D'Val.) V 

if (((hxpiv < 15.0) || (hxtlv < D) || (hxdpv < 1)) 
i 

r = 1.6; 
return(r); 

> 

hxdiv = 1.5260; 
hxdtv = 0.7183; 
hxc = 0.9870; 

betav = hxdtv/hxdiv; /* 0.470707732 */ 

hxrov = (hxpiv - 200.0) / 150.0 * 

(2.096434 - 3.12524E-03*hxt1v + 1 .49281 1E-06*(hxt1v*hxt1v)) ♦ 

(2.728433 - 4.00388E-03*hxt1v + 1 .893453E-06*(hxt1v*hxt1v)); 
hxgamv = (hxpiv - 200.0) / 150.0 * 

(1.007482E-01 - 2.10756E-04*hxt1v + 1 . 14232E-07*(hxt1 v*hxt1v)) ♦ 

(1.543231 - 2.76844E-04*hxt1v + 9.367208E-08*(hxt1v*hxt1v)); 
hxprat = (hxpiv - hxdpv)/hxpiv; 

hxb4 = pow( betav, 4.0); /* 0.049091389 */ 

hxkc = hxc * sqrt(1 .0/(1 .0 - hxb4)); /* 1.01215672 */ 

gfl = (hxgamv - 1.0) /hxgamv; 

gf2 = 2.0/hxgarav; 

sum = pow(hxprat,gf 1); 

gf3 = (1.0 - sun)/(1 .0 - hxprat); 


sum = pow( hxprat ,gf 2); 

gf4 = (1 - hxb4)/(1 - hxb4 * sum); 
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r 1.51897841 V 


hxy * sqrt< sun / gfl * flf3 * gf4); 
hxtxf = 0.99539 ♦ 8.6957E-06 * hxtlv; 
rl = sqrt(2.0 * 32.176 * hxrov * hxdpv); 
r3 « pi * (hxdtv * hxdtv>/4.0/12.0; 
r = hxkc * r3 * hxy * hxtxf * rl; 

return(r); 

} /* woax V 


float wfaxO 

< 

/* Calculate fuel backpressure. Mechanics unknown and assured to work. */ 
/* Robert Van den Heuvel (RSA programmer) 11/26/91 V 

float fpiv, ftlv, fdpv, fc, fbeta, fgam, fgl; 
float fg2, fg3, fg4, fy; 

float fdiv, fdtv, fprat, fb4, fkc, ftxf, frov; 
float pi,r, r1,r2; 
float sun; 

pi * 3.14159; 

fpiv * ACT [53] + 14.7; 
ftlv * ADTK4); 
fdpv = AOT [55] ; 

/* Set UFAX parameters to reasonable values if the measured parameters V 
/* are missing in the incoming data stream <12/5/91 modif. by D'Val): */ 

if (((fpiv < 15.0) || (ftlv < 1)) || (fdpv < 1)) 

< 

r * 0.7; 
return(r); 

y 

fdiv = 1.5200; 
fdtv = 0.7663; 
fc = 0.9870; 

fbeta = fdtv/fdi v; /* 0.504144736 */ 

fgam = 1.365 + 1.6E-4 * (ftlv-400.0) + 1.75E-5 * (fpiv - 400.0); 
fprat = (fpiv - fdpv)/fpiv; 

fb4 * pow( fbeta, 4.0); /* 0.064615194 *f 

fkc = fc * sqrt(1 .0/(1 .0 - fb4)); /* 1.02052112 V 
fgl = (fgam - 1.0)/fgam; 
fg2 = 2.0/fgam; 

sun = pow( fprat, fgl); 


fg3 * (1.0 - sun)/(1 .0 - fprat); 
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sum = pow(fprat,fg2); 


fg4 = <1 - fb4>/(1- fb4 * sun); 

fy = sqrt( sun / fgl * fg3 * fg4); 

ftxf = 0-99539 ♦ 8.6957E-06 * ftlv; /* 1.513006017 */ 

frov = (fpiv - 250.0) / 200.0 * 

(2.457675E-01 • 5.48004E-04*f tlv ♦ 4.07272E-07*(ft1v*f tlv) ) ♦ 
(3.102394E-01 - 6.92109E-04*f tlv ♦ 5.140363E-07*(f t1v*ft1v)); 

rl = sqrt(2.0 * 32.176 * frov * fdpv); /* 15717 */ 
r2 = pi * (fdtv * fdtv)/4. 0/12.0; 

r = fkc * r2 * fy * ftxf * rl; 

return(r); 

> /* wfax V 
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